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This  report  presents  the  results  of  a continued  evaluation  of 
the  Seismic  Research  Observatory  stations  performed  by  Texas  Instruments 
Incorporated  at  the  Seismic  Data  Analysis  Center  in  Alexandria,  Virginia. 
During  the  time  period  this  evaluation  was  carried  out,  data  were  available 
from  six  stations  - Albuquerque,  New  Mexico  (ANMO);  Guam,  Marianas 
Islands  (GUMO);  Mashhad,  Iran  (MAIO),  Narrogin,  Western  Australia  (NWAO); 
Taipei,  Taiwan  (TATO);  and  Wellington  (South  Karori),  New  Zealand  (SNZO). 


The  major  areas  of  investigation  in  the  evaluation  were: 
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• Analysis  of  RMS  noise  levels,  RMS  noise  trends,  and  noise 

spectral  content. 


• Estimation  of  the  quality  of  performance  of  the  short- period 

automatic  event  detector.  * 

• Estimation  of  the  detection  capability  of  the  individual  stations 

on  a regionalized  basis.  ! 

• Estimation  of  the  discrimination  capability  of  the  individual  | 
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stations.  > 

• Estimation  of  the  detection  capability  of  the  stations  functioning 
as  a network. 

In  the  course  of  processing  and  analyzing  data  for  these  areas 
of  investigation,  estimates  were  made  of  the  data  quality,  station  reliability, 
and  probability  of  mixed  events.  Conclusions  regarding  the  above  areas  of 
interest  and  plans  for  further  work  are  also  presented  in  this  report. 
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SECTION  I 
INTRODUCTION 

A.  THE  SEISMIC  RESEARCH  OBSERVATORY  SYSTEM 

It  has  been  noted  (Robinson,  1967)  that  seismic  data  recorded 
by  surface- sited  instruments  is  degraded  or  obscured  to  a significant  degree 
by  wind-induced  earth  tilts.  Theoretical  data  and  tests  indicated  that  this 
wind-induced  noise  component  decreases  rapidly  with  depth.  Therefore,  it 
is  possible  to  enhance  the  seismic  data  quality  by  locating  the  seismic  sensor 
at  a depth  of  approximately  100  meters.  The  Seismic  Research  Observatories 
form  the  first  network  of  single  sites  to  implement  this  observation. 

The  Seismic  Research  Observatories  data  acquisition  and  re- 
cording system  has  been  described  in  detail  previously  (Strauss,  1976). 
Therefore,  only  a brief  description  of  the  system  is  presented  here. 

Broadband  seismic  energy  is  recorded  by  force-balance  type 
seismometers  which  produce  an  output  proportional  to  earth  acceleration  over 
the  frequency  range  of  0.  02  to  1.0  Hz.  Both  long-period  and  short-period  data 
are  produced  from  each  sensor  by  selectively  filtering  the  broadband  output. 
The  long-period  data  are  digitized  and  recorded  continuously  on  an  800  bit  per 
inch  magnetic  tape.  The  vertical  component  of  short- period  motion  is  input 
to  an  automatic  detector,  which  permits  recording  of  vertical-component 
short-period  motion  only  when  certain  operator- specified  conditions  are  sat- 
isfied. Each  magnetic  tape  can  hold  up  to  fourteen  days  of  digital  seismic 
data.  When  full,  these  magnetic  tapes  are  shipped  to  the  Albuquerque  Seis- 
mological  Center  at  Kirtland  Air  Force  Base  in  New  Mexico.  Copies  of  these 
tapes  are  sent  to  the  Seismic  Data  Analysis  Center  at  Alexandria,  Virginia. 


The  long- period  data  are  sampled  once  per  second  and  quantized 
at  5 computer  coxmts  per  millimicron  of  ground  motion.  Short-period  data  are 
sampled  20  times  per  second  and  quantized  at  2000  computer  counts  per  milli- 
micron of  ground  motion  with  the  following  exceptions.  Beginning  1 May  1976 
at  Guam,  14  April  1976  at  Wellington,  New  Zealand,  and  13  May  1976  at  Taipei, 
Taiwan,  the  short-period  data  were  quantized  at  2 computer  counts  per  milli- 
micron of  ground  motion  to  prevent  data  clipping. 

B.  THE  EVALUATION  TASK 

This  report  presents  the  results  of  a continued  evaluation  of  the 
Seismic  Research  Observatories.  The  specific  goals  of  this  evaluation  are: 

e To  estimate  the  data  quality  and  reliability  of  each  station 

• To  investigate  the  short-period  and  long-period  noise  field 
characteristics  of  each  station 

e To  investigate  the  design  of  the  short-period  automatic  detector 

and  to  make  recommendations  to  improve  this  design 

e To  estimate  the  detection  capability  of  each  station  on  a region- 

alized basis 

• To  estimate  the  discrimination  capability  of  each  station  on  a 
regionalized  basis 

• To  investigate  the  performance  of  the  stations  as  a network. 

The  approach  to  reaching  these  evaluation  goals  is  as  follows. 
First,  a suitable  suite  of  seismic  events  and  noise  samples  was  constructed 
from  available  event  lists.  The  procedure  for  selecting  this  event  suite  is 
described  in  Section  II.  This  section  also  describes  the  manner  in  which  the 
events  and  noise  samples  were  processed  and  gives  estimates  derived  from 
processing  experience  of  the  quality  of  the  Seismic  Research  Observatory 
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data  and  reliability  of  the  individual  stations.  Section  III  discusses  the  short- 
period  automatic  detector  resident  at  each  station.  The  discussion  of  this 
topic  covers  the  present  design  of  the  automatic  detector,  design  problems 
found  during  the  course  of  evaluation  of  short-period  data,  and  steps  which 
may  be  taken  to  alleviate  these  problems. 

In  Section  IV,  the  local  noise  field  characteristics  at  each 
evaluated  SRO  station  are  presented.  The  local  noise  field  is  characterized 
by  its  RMS  noise  levels,  monthly  RMS  noise  level  trends,  peak  noise  ampli- 
tudes, and  spectral  content  of  the  noise.  Section  V presents  the  detection 
capability  of  each  station  in  terms  of  the  entire  area  of  interest  and  in  terms 
of  specified  regions  within  the  area  of  interest.  In  the  course  of  estimating 
detection  capability,  the  effect  which  mixed  events  and  system  malfunctions 
have  on  these  estimates  are  discussed.  Section  VI  considers  the  question  of 
discriminating  between  earthquakes  and  presumed  nuclear  explosions.  The 
technique  used  in  this  section  is  the  discriminant.  Section  VII  pre- 

sents a first  look  at  the  detection  capability  of  the  stations  functioning  as  a 
network.  Network  detection  capability  estimates  are  made  from  the  station 
noise  characteristics  and  from  the  detection  statistics  of  the  individual  sta- 
tions. Section  VIII  summarizes  the  results  of  this  evaluation  and  suggests 
future  avenues  of  investigation  to  improve  these  results.  Section  IX  lists 
the  references  cited  in  the  text  of  the  report. 
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SECTION  n 
THE  DATA  BASE 
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A.  DATA  AVAILABILITY 


The  Seismic  Research  Observatory  network,  when  fully  oper- 
ational, will  consist  of  at  least  the  ten  stations  listed  in  Table  II- 1 and  shown 
in  Figure  II- 1.  (The  shaded  area  of  Figure  II- 1 shows  the  area  in  which 
events  used  in  this  evaluation  occurred.  ) During  the  contract  period  covered 
by  this  report,  data  were  available  from  six  of  these  ten  stations:  Albuquerque, 
New  Mexico  (ANMO),  Mashhad,  Iran  (MAIO),  Guam,  Marianas  Islands  (GUMO), 
Narrogin,  Western  Australia  (NWAO),  Taipei,  Taiwan  (TATO),  and  Welling- 
ton (South  Karori),  New  Zealand  (SNZO). 

Figure  II-2  shows  diagrammatically  the  availability  of  data 
for  each  evaluated  station  from  the  start-up  date  of  each  to  the  end  of  1976. 

The  vertical  dashed  lines  in  Figure  II-2  enclose  the  time  frame  from  which 
data  were  selected  for  use  in  performing  this  evaluation.  ANMO,  MAIO,  and 
NWAO  show  near- continuous  data  availability  from  their  start-up  dates  to  the 
end  of  this  time  frame.  GUMO  and  SNZO  each  show  two  major  breaks  in  data 
availability  during  this  time  frame.  Of  the  six  stations,  the  least  data  were 
available  from  TATO,  which  became  operktional  13  May  1976,  was  down  for 
six  days  in  July  1976,  and  went  down  for  the  rest  of  the  year  on  20  September 
1976.  (TATO  did  not  produce  more  data  until  18  July  1977.  ) For  this  reason, 
the  evaluation  time  frame  for  TATO  is  13  May  1976  to  20  September  1976. 

B.  FORMATION  OF  THE  DATA  BASE 

In  this  evaluation  of  the  Seismic  Research  Observatory  stations, 
the  detection  and  discrimination  capabilities  of  each  station  are  to  be  presented 
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on  a regionalized  basis.  The  general  area  of  interest  for  this  evaluation  was 
identified  as  Eurasia.  This  cutback  from  the  world-wide  area  used  in  the 
preliminary  Seismic  Research  Observatory  evaluation  (Strauss,  1976)  was 
made  to  keep  the  size  of  the  data  base  within  reasonable  bounds.  Since 
Eurasia  is  and  always  has  been  the  main  area  of  interest  from  an  evaluation 
standpoint,  this  cutback  can  be  made  without  any  significant  loss  of  information. 

The  first  step  in  forming  the  evaluation  data  base  was  to  se- 
lect all  seismic  events  with  Eurasian  epicenters  as  listed  by  the  Norwegian 
Seismic  Array  (NORSAR)  bulletins.  (At  the  time  of  formation  of  the  data 
base,  this  was  the  only  available  event  list.  ) The  time  frame  selected  was 
from  22  December  1975,  when  the  first  station  became  operational,  to  30 
September  1976,  when  the  NORSAR  bulletins  ceased  being  issued.  This 
selection  yielded  a Eurasian  event  list  of  2697  events.  Since  this  was  far 
too  great  a volume  of  data  to  be  processed  and  analyzed  in  the  available  time, 
it  was  necessary  to  reduce  this  list  to  manageable  proportions. 

The  event  list  was  reduced  in  the  following  manner.  First, 
the  events  were  grouped  by  seismic  region  as  defined  by  Flinn  and  Engdahl 
(Flinn  and  Engdahl,  1965).  Those  seismic  regions  containing  fewer  than 
80  events  were  immediately  rejected  and  their  contained  events  dropped  from 
the  event  list.  Next,  the  events  in  the  remaining  regions  were  grouped  by 
their  sub-regions,  again  as  defined  by  Flinn  and  Engdahl.  Those  sub-regions 
which  contained  only  a few  events  and  which  were  separated  from  the  main 
body  of  events  in  that  region  were  then  rejected  and  their  contained  events 
deleted  from  the  event  list. 

Table  U-2  lists  the  region  numbers  and  corresponding  sub- 
region  numbers  for  the  data  base.  Flinn  and  Engdahl' s region  19  was  sub- 
divided into  regions  1 and  2,  corresponding  to  the  Kurile  Islands  and  northern 
Japan,  respectively.  This  sub-division  was  made  to  centralize  the  grouping 
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of  events.  The  regions  were  re-numbered  simply  to  make  them  sequential 
and  will  be  referred  to  by  these  new  numbers  throughout  this  report.  (Not 
all  of  these  regions  are  used  in  the  discussion  of  regionalized  results  of 
Sections  V and  VI  since  mixed  events  and  system  malfunctions  in  some  cases 
reduced  the  number  of  events  in  a region  below  the  point  where  reasonable 
results  could  be  determined.  ) The  region  boundaries  are  shown  in  Figure 
n-3. 
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At  this  point,  the  data  base  still  contained  approximately  2300 
events.  To  further  reduce  this  number,  events  were  systematically  deleted 
(every  other  event,  two  out  of  three,  as  needed)  from  heavily  event-populated 
regions  until  no  region  contained  more  than  150  events.  This  procedure  re- 
sulted in  the  final  data  base  of  1472  events.  The  event  parameters  of  the 
data  base  are  given  in  Appendix  A. 

It  should  be  noted  that  the  manner  in  which  the  data  base  was 
created  produces  a data  base  which  is  essentially  a random  sampling  of  Eur- 
asian seismic  events.  No  effort  was  made  to  exclude  obviously  mixed  events, 
so  that  mixed  event  statistics  derived  from  this  data  base  should  accurately 
reflect  the  frequency  of  occurrence  of  mixed  events.  The  frequency  of  occur- 
ence of  system  malfunctions  resulting  in  the  inability  of  the  analyst  to  study 
a desired  segment  of  seismic  data  will  yield  representative  estimates  of  sta- 
tion reliability. 

The  short-period  data  base  used  in  this  evaluation  is  simply 
a subset  of  the  above -described  data  base.  All  members  of  this  subset  have 
station- epicenter  separations  of  less  than  103  degrees.  This  restriction  was 
placed  on  the  short-period  data  base  since  detection  of  the  P wave  beyond 
103**  from  the  event  epicenter  is  limited  by  a shadow  zone.  At  some  stations 
this  restriction  severely  limited  the  size  of  the  data  base,  since  these  sta- 
tions are  located  far  from  the  general  area  of  interest.  Table  II- 3 presents 
the  mean  epicentral  distances  from  each  station  to  each  region.  (These  values 
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are  derived  from  the  station- epicenter  separations  of  the  events  occurring  in 
each  region.  ) In  particular,  note  that  station  SNZO  is  more  than  103  degrees 
from  the  majority  of  the  regions.  Therefore,  very  little  can  be  expected  from 
the  short-period  data  recorded  at  this  station. 

The  description  of  the  noise  sample  data  base  is  given  in  detail 
in  Section  IV  of  this  report.  In  brief,  the  short-period  noise  samples  were 
selected  from  time  gates  immediately  preceding  the  observed  signal  which 
triggered  the  automatic  detector.  Long-ptriod  noise  samples  were  selected 
from  time  gates  during  which  no  signals  were  expected  to  be  arriving  at  the 
stations  as  determined  from  the  NORSAR  event  bulletins. 

C.  DATA  PROCESSING 

1.  Long-Period  Data  Processing 

The  long-period  signal  and  noise  data  used  in  this  evaluation 
were  processed  in  two  stages--a  basic  processing  stage  and  an  analysis  pro- 
cessing stage.  The  basic  processing  stage  is  illustrated  by  Figure  II-4.  For 
long-period  data,  the  basic  processing  stage  exists  in  two  forms.  In  the  first 
form,  the  seismic  events  and  noise  samples  selected  for  analysis  are  first 
processed  by  an  edit  program,  which  searches  the  field  tape  for  the  desired 
edit  start  time,  reads  in  the  desired  data,  resamples  the  data  to  a 2- second 
sample  rate,  quality  checks  the  data,  computes  the  trace  means,  and  writes 
the  data  on  an  edit  tape.  At  this  point,  the  three  components  of  motion  are 
in  their  recorded  vertical,  north,  east  (V,  N,  E)  configuration.  Next,  the 
data  is  read  from  the  edit  tape  by  a program  which  removes  the  trace  means, 
rotates  the  data  from  its  V,  N,  E configuration  to  a vertical,  transverse, 
radial  (V,  T,  R)  configuration,  and  writes  the  result  on  a second  tape. 

Finally,  this  rotated  data  is  read  by  a third  program  which  bandpass  filters 
(0.  023-0.  059  He  passband)  and  plots  the  data. 
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This  processing  scheme  was  originally  designed  for  the  evalu- 
ation of  the  high-gain  long-period  (HGLP)  stations,  also  known  as  the  Very 
Long -Period  Experiment  (VLPE)  stations.  With  the  inception  of  the  evaluation 
of  the  Seismic  Research  Observatories,  the  VLPE  programs  were  modified 
to  handle  the  new  data  formats.  When  these  programs  were  originally  written, 
it  was  felt  that  a great  deal  of  analyst  interaction  with  the  data  would  be  neces- 
sary. Thus,  the  processing  package  was  deliberately  designed  in  the  above- 
described  separate  levels  of  edit,  rotate,  end  filter.  These  programs  also 
included  a number  of  options  designed  to  provide  analysis  tools,  such  as 
matched  filters,  which  subsequently  turned  out  to  be  rarely  used.  The  net 
result  was  a fairly  cumbersome  processing  method. 


During  the  current  contract  period,  it  became  obvious  that 
this  form  of  the  basic  processing  stage  was  too  expensive  in  terms  of  analyst 
and  computer  time  to  remain  in  use.  A secondary  problem  with  this  form  of 
the  basic  processing  stage  is  that  it  requires  both  edit  and  rotated  data  hold 
tapes,  when  only  the  rotated  data  tapes  are  necessary,  thus  doubling  magnetic 
tape  consumption.  Therefore,  the  edit,  rotate,  and  filter  programs  were 
combined  as  illustrated  in  Figure  11-4  by  the  area  enclosed  by  the  dashed  line. 
The  basic  processing  stage  in  this  new  form  uses  approximately  one-third  to 
one-half  as  much  computer  time  as  the  old  form,  greatly  reduces  the  time 
which  the  analyst  must  spend  processing  the  data,  and  reduces  tape  consump- 
tion by  one -half. 

Naturally,  some  sacrifices  had  to  be  made  to  build  a program, 
encompassing  the  functions  of  three  programs,  which  would  fit  the  core  lim- 
itations of  the  computer.  Therefore,  the  least-used  options  of  the  old  package 
were  dropped  when  the  new  program  was  formed.  These  options  have  not 
been  lost,  however.  Should  the  analyst  find  the  need  for  any  of  them,  he  can 
revert  to  the  use  of  the  original  programs,  which  have  been  retained  in  the 
form  of  Fortran  card  decks  and  listings. 
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The  new  program  was  thoroughly  checked  before  it  replaced 
the  old  programs.  The  method  of  checking  was  to  perform  a hexadecimal 
dump  of  a portion  of  a fteld  tape.  These  data  were  manually  decoded  from 
hexadecimal  to  decimal  and  the  appropriate  trace  means  subtracted.  The 
data  were  next  manually  rotated  from  their  V,  N,  E configuration  to  the 
V,  T,  R configuration.  Finally,  the  results  were  compared  to  the  equivalent 
data  output  from  the  new  program.  The  two  sets  of  data  were  found  to  agree 
in  all  particulars.  As  a final  check,  several  events  were  processed  using  the 
old  programs  and  the  new  program.  The  results  were  compared  and  found 
to  be  identical. 

One  other  benefit  of  the  new  basic  processing  stage  was  found 
after  it  was  put  in  use.  In  the  old  edit  program,  if  a parity  error  was  encoun- 
tered while  the  program  was  searching  for  the  edit  start  time,  the  program 
would  be  unable  to  proceed  and  the  analyst  was  forced  to  declare  a malfunction 
(tape  error)  for  that  event.  To  avoid  this  problem,  when  the  new  program 
encounters  a parity  error,  it  increments  a parity  counter  and  proceeds.  The 
total  number  of  parity  errors  is  printed  out  after  the  edit  start  time  is  found. 
This  approach  has  greatly  reduced  the  number  of  declared  malfunctions.  Doc- 
umentation for  this  new  basic  processing  stage  has  been  completed  during  the 
current  contract  period  (Schmidt,  1977). 

The  long-period  analysis  processing  stage  is  described  in  de- 
tail in  Section  IV  for  noise  analysis,  in  Section  V for  signal  detection,  and  in 
Section  VI  for  signal  discrimination.  Therefore,  this  stage  will  be  only  briefly 
described  here. 

Once  the  long-period  noise  samples  have  been  processed  through 
the  basic  processing  stage,  they  are  visually  checked  for  unreported  events 
and  system  malfunctions  such  as  spikes  or  glitches.  Each  noise  sample  pas- 
sing this  visual  check  is  then  processed  by  the  analysis  program,  which  com- 
putes RMS  noise  amplitudes,  measures  the  peak  25- second  noise  amplitude  of 
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the  sample,  and  computes  the  power  spectrum  of  each  component  of  the  sample. 
All  computed  values  are  output  on  punched  cards,  which  form  the  input  to  the 
programs  which  plot  RMS  noise  amplitudes,  RMS  noise  trends,  spectral  con- 
tent of  the  noise,  and  maximum  25- second  noise  amplitudes. 

The  signal  analysis  stage  is  primarily  visual,  as  the  analyst 

searches  each  record  for  detection  of  the  event  under  consideration  while 

simultaneously  checking  for  data  quality.  One  program  used  by  the  analyst 

in  this  stage  is  the  maximum  likelihood  curve  fitting  program  (Ringdal,  1974) 

which  is  used  in  estimating  station  detection  capability.  A second  program 

used  in  the  signal  analysis  is  the  surface  wave  magnitude  (M^)  measurement 

program.  This  program  computes  the  surface  wave  magnitude  for  each 

waveform  in  the  signal  gate  and  then  searches  this  series  of  M values  for 

8 

the  largest  20,  25,  and  30  second  M values. 

8 

2.  Short-Period  Data  Processing 

Unlike  the  long-period  data  processing,  the  short-period 
data  processing  was  performed  in  a single  stage.  This  processing  was 
carried  out  on  the  PDP-15  computer  using  the  Interactive  Seismic  Processing 
System  (ISPS)  package  (Ringdal,  Shaub,  and  Black,  1975).  The  processing 
method  is  shown  in  Figure  II- 5. 

The  short-period  record  headers  on  the  1600  BPI  copy  of  the 
SRO  field  tape  are  first  dumped  to  reveal  when  the  automatic  detector  permit- 
ted the  recording  of  short-period  data.  The  data  base  is  then  searched  to  find 
those  events  for  which  short- period  data  were  recorded.  These  events  are  then 
edited,  resampled  to  10  points  per  second,  and  written  on  the  disk  pack.  The 
analyst  can  then  filter  the  data,  (0.  5-4.  0 Hz  passband)  make  Calcomp  plots, 
and  measure  such  parameters  as  RMS  noise  amplitudes  and  peak  one- second 
noise  amplitudes.  When  the  analyst  is  finished  with  the  data  from  one  tape, 
the  disk  data  files  are  overwritten  with  data  from  the  next  tape  and  analysis 
is  resumed. 
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FIGURE  n-5 

SHORT-PERIOD  DATA  PROCESSING  FLOW  CHART 
(PDF- 15  DATA  PROCESSING) 
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Obviously,  the  major  drawback  to  this  data  processing  approach 
is  that  the  processed  short-period  data  is  not  saved  on  a permanent  hold  tape. 
Thus,  if  the  analyst  finds  it  necessary  to  re-examine  a particular  data  seg- 
ment, he  must  repeat  the  entire  processing  sequence.  A.  second  drawback 
to  this  method  of  short-period  data  processing  was  the  unreliability  of  the 
peripheral  units  of  the  PDP-15.  Malfunctioning  tape  drives  stopped  short- 
period  data  processing  completely  for  two  months.  Other  peripheral  units 
which  were  subject  to  failure  (albeit,  to  a lesser  degree)  were  the  line  printer 
and  the  plotter. 

In  order  to  provide  a permanent  short-period  data  hold  capabil- 
ity and  to  provide  a back-up  processing  system  for  those  times  when  the  PDP- 
15  is  down,  a short-period  edit  program  for  use  on  the  IBM  360/44  was  created 
late  in  the  contract  period  (Schmidt,  1977).  The  edited  data  are  written  on  tape 
in  the  same  format  as  that  of  the  long-period  rotated  data.  A program  has 
also  been  created  to  filter  and  plot  the  data  and  measure  desired  parameters. 


are: 


The  products  of  the  short-period  and  long-period  processing 

Permanent  hold  tapes  containing  processed  signal  and  noise 
data 

Calcomp  plots  of  the  bandpass  filtered  data 
RMS  noise  measurements 
Noise  spectra 

Maximum  1 second  (SP)  and  25  second  (LP)  noise  amplitudes 
Detection  statistics 

Surface  wave  magnitudes  computed  at  20,  25,  and  30  second 
periods. 
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I. 

I D.  PROCESSING  SUMMARY 

Table  II-4  summarizes  the  processing  results  of  the  work  per- 
I formed  during  the  current  contract  period  on  the  Seismic  Research  Observa- 

tories evaluation  task.  In  this  table,  the  'SP'  and  'LP'  under  the  heading 
I 'DATA  TYPE'  refer  to  short-period  data  and  long-period  data  respectively. 

The  heading  'DETECTED'  refers  to  the  n\imber  of  events  which  were  visually 
j detected  on  a bandpass  filtered  (0.  023-0.  C59  Hz  passband  for  LP  and  0.  5-4.  0 

Hz  passband  for  SP)plot  under  the  detection  criteria  of  Section  V.  The  heading 
{ 'NOT  DETECTED'  refers  to  the  number  of  events  for  which  only  seismic  noise 

was  observed  in  the  signal  gate.  The  heading  'MIXED'  refers  to  the  number  of 
( events  which  were  not  detectable  due  to  the  presence  of  some  other  signal  in  the 

signal  gate.  The  heading  'NO  SP  DATA  RECORDED'  refers  to  the  number  of 
I events  not  detected  by  the  short-period  automatic  detector.  The  heading  'SYS- 

TEM FAILURE-NO  DATA  RECORDED'  refers  to  the  number  of  events  for  which 
( no  data  were  recorded  since  the  system  had  gone  down  (failed).  The  heading 


'MALFUNCTION'  refers  to  the  number  of  events  for  which  the  detection  sta- 
tus could  not  be  determined  due  to  malfunctions  (spikes,  glitches,  and  the 
like)  in  the  signal  gate.  The  heading  l^<103*^'  refers  to  the  number  of  events 
for  which  short-period  data  were  not  processed  due  to  their  epicentral  dis- 
tances being  greater  than  103°. 

It  was  not  possible  to  process  all  of  the  events  in  the  data 
base  at  all  of  the  stations  under  evaluation.  Conversion  of  the  long-period 


processing  package  and  the  subsequent  checking  out  of  the  new  processing 

I 

I I program  made  substantial  inroads  into  the  time  available  for  data  processing. 

Short- period  data  processing  was  delayed  by  hardware  failures  on  the  PDP-15. 
I I Even  with  these  difficulties,  however,  a substantial  amount  of  data  was  pro- 

i cessed  and  analysed  during  the  current  contract  period. 
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The  quality  of  the  data  recorded  at  the  Sesimic  Research  Ob- 
servatories was  very  good  in  most  cases.  The  types  of  malfunctions  which 
degrade  the  data  quality  were  described  in  detail  in  the  preliminary  evaluation 
of  the  Seismic  Research  Observatories  (Strauss,  1976).  The  overall  frequency 
of  occurrence  of  these  malfunctions  can  be  estimated  from  Table  II-4.  Long- 
period  data  were  degraded  by  malfunctions  three  to  eight  percent  of  the  time. 
With  the  exception  of  data  recorded  at  MAI.O,  the  short-period  data  is  degraded 
by  malfunctions  one  percent  of  the  time  or  less.  The  relatively  high  malfunc- 
tion rate  of  MAIO  short-period  data  is  primarily  due  to  a type  of  malfunction 
first  encountered  during  the  current  contract  period.  This  malfunction  appears 
on  the  short-period  record  as  a 1.  7 Hz  (approximate)  waveform.  The  waveform 
rises  out  of  the  noise  field  with  an  onset  time  of  approximately  three  seconds. 
This  waveform  reaches  peak  amplitudes  (measured  zero-to-peak)  of  four  to 
five  millimicrons  and  may  continue  for  several  hours.  (The  effect  of  this 
on  the  short-period  automatic  detector  is  discussed  in  Section  III.  ) 

Estimates  of  station  reliability  can  be  made  from  the  values 
given  in  Table  II-4.  The  two  factors  which  detract  from  the  reliability  of  the 
station  are  the  number  of  times  it  was  not  possible  to  examine  the  data  since 
no  data  was  recorded  and  the  number  of  times  it  was  not  possible  to  examine 
the  data  due  to  system  malfunctions.  Thus,  the  reliability  estimates  presented 
here  are  calculated  from: 

Reliability  = 1 - ^(no  data  recorded)  + P(malfunction)j 

where  the  two  probabilities  are  estimated  empirically  from  the  data  of  Table 
n-4. 

The  reliability  estimates  are  presented  in  Table  II-5.  The  rel- 
atively low  reliability  of  GUMO  is  primarily  due  to  station  down  time,  when 
no  data  was  recorded  (see  Figure  11-2). 
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TABLE  n-5 

STATION  RELIABILITY  ESTIMATES 


Station 


ANMO 

GUMC 

MAIO 

NWAO 

TATO 

SNZO 


Reliability 

Short-Period 

Long-Period 

C.  91 

0.87 

0.76 

0.70 

0.89 

0.  92 

1.00 

0.91 

0.  94 

0.  90 

0.94 

0.80 

2 


• Stable  false  alarm  rate  over  a wide  range  of  ambient  noise 
levels 

e Minimum  number  of  detections  declared  from  the  same  event 

while  retaining  the  ability  to  detect  possible  secondary  phase 
arrivals  within  the  coda 

e Minimum  time  necessary  to  re-establish  equilibrium  after 

parameter  value  changes 

• Low  sensitivity  to  degraded  data  quality  (Swindell  and  Snell, 
1977). 

In  this  evaluation  of  the  automatic  detector  design,  data  output 
by  the  automatic  detector  will  be  assessed  in  terms  of  the  above  points. 


four  events 


e 


e 


At  any  point  in  the  operation  of  the  automatic  detector,  one  of 
must  occur: 

Signal  in  time  gate,  detection  declared 

No  signal  in  time  gate,  no  detection  declared 

Signal  in  time  gate,  no  detection  declared  (type  1 error) 

No  signal  in  time  gate,  detection  declared  (type  2 error). 


The  third  and  fourth  of  these  four  events  are,  of  course,  the  two  possible 
types  of  detector  errors.  With  only  data  which  has  been  processed  by  the 
automatic  detector  available,  it  is  not  possible  to  assess  the  frequency  of 
occurrence  of  the  type  1 errors.  It  would  be  necessary  to  have  continuous 
data  and  data  for  the  same  time  frame  which  has  been  processed  by  the  auto- 
matic detector  to  determine  how  often  the  automatic  detector  fails  to  detect 
an  analyst-detectable  signal.  It  is  possible  to  at  least  qualitatively  assess 
the  performance  of  the  automatic  detector  in  terms  of  the  type  2 error,  i.  e. , 
false  alarms. 
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B.  DESIGN  OF  THE  PRESENT  SHORT- PERIOD  AUTOMATIC  DETECTOR 

Each  Seismic  Research  Observatory  station  is  equipped  with  a 
short-period  automatic  event  detector  (hereafter  referred  to  as  the  SRO  auto- 
matic detector)  which  is  intended  to  discriminate  between  seismic  noise  and 
short-period  seismic  events.  The  design  of  the  SRO  automatic  detector  is 
based  on  a study  carried  out  by  the  Charles  Stark  Draper  Laboratory  (Eterno, 
etal. , 1974). 

The  design  concept  of  the  SRO  automatic  detector  is  illustrated 
by  Figure  IIl-l.  The  short- period  analog  data  from  the  sensor  unit  is  first 
digitized  at  a rate  of  twenty  points  per  second.  At  the  SRO  automatic  detector 
input,  each  datum  is  converted  to  a twelve  bit  integer  quantity.  The  digitized 
data  next  passes  through  a band-limiting  four  pole-pair  Butterworth  filter 
with  two  zeros  at  zero  frequency.  (This  filter  is  constructed  by  cascading 
four  recursive  single  pole- pair  filters.  A single- zero  filter  is  placed  in  front 
of  the  first  and  fourth  cascaded  filters.  ) The  lower  cutoff  of  this  filter  is  0.  5 
Hz,  the  upper  cutoff  is  2.  0 Hz,  the  filter  response  within  the  passband  is  flat 
to  within  ^ 1 dB,  and  the  rejection  rate  outside  the  passband  is  more  than  25 
dB. 

The  data  are  then  rectified  and  passed  to  the  short-term  aver- 
ager (STA).  Utilizing  a single-pole  recursive  stage,  one- second  exponentially - 
weighted  short-term  averages  are  constantly  generated.  The  short-term  aver- 
ager imparts  an  absolute  gain  of  approximately  13  dB  to  the  data.  The  short- 
term averages  are  combined  to  form  a long-term  average  (LTA).  Long-term 
averages  may  be  formed  over  operator- selected  time  intet-'als  of  10,  20,  30, 
or  40  seconds.  The  approximate  13  dB  gain  imparted  by  the  short-term  aver- 
ager allows  greater  accuracy  in  calculating  the  long-term  average.  However, 
once  the  calculation  is  complete,  the  long-term  average  is  rescaled  to  unitize 
the  gain  at  zero  frequency. 
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c. 


PROBLEMS  WITH  THE  CURRENT  DESIGN  OF  THE  SRO  AUTOMATIC 
DETECTOR 


During  the  course  of  processing  and  analyzing  short>period  data 
to  be  used  in  the  estimation  of  the  short-period  detection  capability  of  the  Seis- 
mic Research  Observatory  stations,  several  persistent  problems  were  found 
in  the  present  design  of  the  SRO  automatic  detector. 

The  first  problem  found  in  the  automatic  detector  design  is 
illustrated  by  Figure  1II-2.  The  monotonic  waveform  (1.7  Hz  frequency  ap- 
proxinnately)  seen  emerging  from  the  noise  in  the  upper  trace  has  been  re- 
peatedly found  to  be  the  triggering  agent  of  the  SRO  automatic  detector  at 
ANMO  and  MAIO.  The  source  of  this  waveform  at  ANMO  is  apparently  a 
large  centrifuge  operated  at  Kirtland  Air  Force  Base,  while  at  MAIO  it  is 
apparently  a rock  crusher.  The  center  trace  of  Figure  III-2  shows  the  data 
as  it  would  appear  to  the  SRO  automatic  detector  after  application  of  the  0.  5- 
2.  0 Hz  passband  filter.  This  filter  has  cleaned  up  the  seismic  noise  preceed- 
ing  the  1.  7 Hz  filter  while  leaving  this  waveform  essentially  xineffected.  The 
SRO  automatic  detector  detects  the  rise  in  the  short-term  average  due  to 
the  onset  of  this  waveform  and  starts  declaring  detections.  Since  this  wave- 
form at  times  persists  for  hours,  hundreds  of  false  alarms  can  be  rapidly 
written  on  tape. 

The  lower  trace  of  Figure  HI- 2 presents  the  same  data  after 


application  of  a 0.  5-1.  5 Hz  passband  filter.  Note  that  the  1.  7 Hz  waveform 
has  been  completely  removed  from  the  data.  The  scale  factor  for  this  trace 
indicates  that  the  remaining  waveforms  are  simply  seismic  noise.  Tests  of 
this  filter  on  seismic  signals  buried  in  this  1.7  Hz  waveform  indicate  that  i 


the  1. 7 Hz  waveform  can  be  removed  without  significantly  degrading  the  sig- 
nal. Therefore,  the  detection  false  alarm  rate  can  be  substantially  reduced 
at  stations  displaying  this  1.  7 Hz  waveform  by  changing  the  passband  filter 
limits  of  the  SRO  automatic  detector  from  0.  5-2.  0 Hz  to  0.  5-1. 5 Hz. 
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FIGURE  m-2 

SHORT- PERIOD  NOISE  SAMPLE  SHOWING 
ONSET  OF  1 . 7 Hz  WAVEFORM 
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The  next  design  problem  to  be  discussed  is  that  the  false  alarm 
rate  of  the  SRO  automatic  detector  is  not  stable.  Frequently,  once  the  detec- 
tor declares  a detection,  it  will  continuously  declare  detections  for  hours.  In 
the  worst  case,  at  SNZO,  short-period  data  has  been  continuously  recorded 
for  more  than  a day  at  a time.  However,  when  the  short-period  data  following 
the  first  short-period  record  of  a continuous  series  are  examined,  in  many 
cases  no  waveform  can  be  found  which  covild  have  caused  the  detector  to  trig- 
ger. From  the  standpoint  of  the  analyst,  these  detections  must  be  called  false 
alarms. 

The  cause  of  these  multiple  detection  declarations  appears  to 
be  involved  with  the  freezing  of  the  long-term  average  when  a detection  is  de- 
clared. This,  coupled  with  the  turn-off  threshold,  causes  the  SRO  automatic 
detector  to  continue  to  declare  detections  until  the  short-term  average  drops 
sufficiently.  The  intent  of  this  arrangement  was  to  permit  the  recording  of 
the  P-wave  coda  and  secondary  body  phases,  hi  practice,  the  result  is  quite 
different.  For  example,  when  the  previously  described  1. 7 Hz  waveform  trig- 
gers the  detector,  the  frozen  long-term  average  and  turn-off  threshold  ensure 
that  short-period  data  recording  will  be  continuous  from  that  point  until  this 
waveform  drops  back  into  the  noise. 

If  the  short-period  noise  level  was  to  abruptly  rise  and  hold  its 
new  level,  a long  series  of  detections  would  be  made.  This  could  happen  if 
the  tum-on  threshold  was  set  too  low  and  the  station  was  sited  near  a coast- 
line, where  an  approaching  storm  could  provide  the  increase  in  noise  level. 
This  may  be  the  cause  of  the  long  series  of  continuously  declared  detections  at 
SNZO,  where  visual  examination  of  the  data  has  revealed  no  reason  for  the 
declared  detections. 

Ehcamples  of  the  instability  of  the  SRO  automatic  detector  can 
be  seen  in  Tables  III-l  and  HI-Z.  Table  ni-l  presents  the  detection  statistics 
for  a 10-day  period  in  1976.  In  this  period,  the  automatic  detector  at  GUMO 
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TABLE  ni-1 

PERFORMANCE  OF  AUTOMATIC  DETECTOR  FOR  AN  ARBITRARY 

10- DAY  PERIOD  IN  1976 


Number  of  Detections 


ANMO  CUMO  MAIO  NWAO  TATO  I SNZO 


06/21/76 

06/22/76 

06/23/76 

06/24/76 

06/25/76 

06/26/76 

06/27/76 

06/28/76 

06/29/76 

06/30/76 


06/21/76 

06/22/76 

06/23/76 

06/24/76 

06/25/76 

06/26/76 

06/27/76 

06/28/76 

06/29/76 

06/30/76 


Percentage  of  Time  Recorded 


ANMO  GUMO  MAIO  NWAO  TATO 


SNZO 
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TABLE  ni-2 

PERFORMANCE  OF  AUTOMATIC  DETECTOR  FOR  AN  ARBITRARY 

9-DAY  PERIOD  IN  1977 

4 


reported  many  more  detections  than  did  the  other  stations.  This,  coupled 
with  the  high  short-period  noise  levels  recorded  at  GUMO  (Section  IV)  leads 
the  authors  to  believe  that  many  of  these  reported  detections  are  false  alarms. 
This  table  also  shows  that  the  percentage  of  short-period  data  recorded  (where 
100  percent  indicates  continuous  recording)  varies  widely  from  day  to  day  and 
from  station  to  station.  For  example,  compare  days  6/21/76  and  6/27/76  at 
ANMO.  Although  the  number  of  detections  declared  are  about  the  same  for 
the  two  days,  the  percentage  of  recorded  data  differs  greatly  (59  percent  as 
opposed  to  2 percent). 

Since  software  and  parameter  setting  changes  have  been  made 
following  this  time  period.  Table  HI- 2 was  constructed  to  illustrate  the  auto- 
matic detector  performance  approximately  one  year  later.  The  data  in  this 
table  indicate  greater  stability  in  detector  performance.  However,  several 
points  should  be  noted.  On  7/29/77  ANMO  shows  67  detections,  a value  far 
greater  than  thos^  at  the  other  stations.  Also,  at  GUMO  on  7/31/77,  97  per- 
cent of  the  short-period  data  were  recorded,  an  amount  that  appears  far  too 
large  in  view  of  the  amounts  of  short-period  data  recorded  at  the  other  stations 
on  that  day.  Finally,  the  lack  of  short-period  data  at  SNZO  indicates  that  per- 
haps changes  in  parameter  settings  went  to  the  opposite  extreme  from  those 
of  a year  previously.  It  would  appear  that  the  detection  thresholds  were  set 
too  high  to  permit  any  detections. 

Next  to  be  described  are  two  sources  of  false  alarms  which, 
while  less  common  than  those  previously  described,  must  be  recognized  and 
corrected.  Figure  III- 3 shows  in  the  upper  trace  an  unfiltered  short-period 
record  from  SNZO  and  in  the  lower  trace  the  same  record  after  filtering  by 
a 0.  5 to  2.  0 Hz  filter.  This  record  was  the  first  in  a series  of  continuous 
detection  declarations.  Therefore,  the  reason  this  record  was  declared  a 
detection  by  the  SRO  automatic  detector  must  be  the  waveform  produced  by 
filtering  the  spike.  This  type  of  false  alarm  can  be  eliminated  by  placing  a 
spike-correction  algorithm  at  the  front  end  of  the  SRO  automatic  detector. 
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FIGURE  m-3 

SHORT-PERIOD  NOISE  SAMPLE  CONTAINING  SPIKE 
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Figure  III-4  shows  the  second  of  these  false  alarm  sources. 
The  waveform  causing  the  SRO  automatic  detector  to  declare  a detection  is 
a glitch  of  unknown  source;  it  does  not  resemble  any  of  the  recorded  wave- 
forms which  the  author  has  been  able  to  associate  with  known  seismic  events. 
In  particular,  it  displays  no  recognizable  coda  following  the  initial  pulse. 

This  waveform  triggered  the  automatic  detector  by  exceeding  the  turn-on 
threshold  for  more  than  twenty  sample  points.  However,  if  it  had  been  re- 
quired that  the  short-term  average  exceed  the  scaled  long-term  average  for 
more  than  twenty  points  (say,  thirty  points)  this  false  alarm  would  not  have 
been  declared. 

In  summation,  the  following  changes  in  the  design  of  the  SRO 
automatic  detector  should  be  made: 

• Change  the  bandpass  filter  limits  from  0.  5-2.  0 Hz  to  0.  5-1. 5 
Hz,  or  notch  out  a very  narrow  band  about  1. 7 Hz. 

• Alter  the  frozen  long-term  turn-off  threshold  combination  so 
that  fewer  detection  declarations  following  the  first  of  a con- 
tinuous series  are  made. 

• Put  a spike -correction  algorithm  at  the  front  end  of  the  SRO 
automatic  detector. 

• Increase  the  number  of  contiguous  points  for  which  the  short- 
term average  must  exceed  the  scaled  long-term  average  in 
order  for  a detection  to  be  declared  from  20  points  to  at  least 
30  points. 
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FIGURE  m-4 

SHORT-PERIOD  NOISE  SAMPLE  CONTAINING  GLITCH 


SECTION  IV 
NOISE  ANALYSIS 


A.  DISCUSSION 

The  goal  of  this  section  is  to  describe  the  characteristics  of 
the  noise  field  at  each  Seismic  Research  Observatory  under  evaluation.  The 
preliminary  Seismic  Research  Observatory  evaluation  (Strauss,  1976)  present- 
ed a first  look  at  short-period  and  long-period  RMS  noise  levels  and  at  long- 
period  average  noise  spectra.  The  long-period  noise  was  corrected  for  in- 
strument response  while  Hie  short-period  results  were  not,  since  instrument 
response  correction  capability  is  not  part  of  the  PDP-15  interactive  processing 
system. 

In  this  continviation  of  the  Seismic  Research  Observatories 
evaluation,  one  change  in  and  one  addition  to  the  noise  analysis  have  been 
made.  First,  both  short-period  and  long-period  RMS  noise  and  spectra  con- 
tent results  are  presented  without  instrument  response  corrections.  This 
change  was  made  since  it  is  recognized  that  the  analyst  is  primarily  concern- 
ed with  the  noise  as  he  will  see  it,  i,  e.,  after  it  has  passed  through  the  sensing, 
filtering,  and  recording  instrumentation.  A long-period  noise  analysis  with 
instrument  response  corrections  is  presented  in  Appendix  B for  the  interested 
reader. 

The  addition  to  the  noise  analysis  is  the  measurement  of  max- 
imum noise  amplitudes,  which  are  made  at  l-second  periods  for  short-period 
data  and  25- second  periods  for  long-period  data.  These  measurements  were 
made  at  the  suggestion  of  Dr.  Filson  of  the  Advanced  Research  Projects 
Agency.  They  will  be  used  in  the  analysis  of  Seismic  Research  Observatories 
network  detection  capability  (Section  VII). 
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The  definition  of  eeismic  noise  used  in  this  report  is  most 
clearly  stated  by  Enders  A.  Robinson  (Robinson,  1967): 

Any  ground  motion  that  is  not  caused  by  an  explosion 
or  an  earthquake  is  usually  regarded  as  ambient  seis- 
mic noise.  The  predominant  components  of  such  seis- 
mic noise  are  surface- generated  microseisms  that 
originate  from  meteorological,  hydrodynamic,  or  cul- 
tural sources.  Such  microseisms  chiefly  propagate 
along  the  surface  of  the  earth  as  Rayleigh  waves. 

B.  VERTICAL-COMPONENT  SHORT-PERIOD  NOISE 

Short- period  noise  analysis  was  performed  on  the  PDP-15 
computer  using  the  Interactive  Seismic  Processing  System  (Ringdal,  Shaub, 
and  Black,  1975).  Subject  to  data  availability  and  time  constraints,  a min- 
imum of  two  noise  samples  per  day  from  each  station  were  processed.  These 
noise  samples  were  obtained  in  one  of  two  ways.  First,  while  processing 
short-period  events  of  the  data  base,  a set  of  noise  samples  was  obtained  by 
selecting  the  time  gate  prior  to  the  observed  P-wave  arrival.  A second  set 
of  short-period  edits  was  made  to  produce  noise  samples  for  days  on  which 
the  data  base  did  not  have  any  events.  Each  of  these  edits  started  at  the  first 
point  of  a data  segment  output  by  the  automatic  detector.  This  start  time 
was  picked  to  avoid  making  noise  measurements  in  signal  coda.  Due  to  the 
manner  in  which  the  automatic  detector  operates,  the  length  of  the  time  gate 
prior  to  the  P-wave  arrival  (and,  hence,  the  length  of  the  noise  sample)  was 
usually  about  20  seconds. 

After  the  noise  gate  was  selected,  each  sample  was  filtered  by 
a 0.  5-4.  0 Hz  bandpass  filter.  RMS  noise  was  then  computed  by  the  equation 

tl  '**’"1  * 

RMS  Noise  » ^ 


IV- 2 


► 


I 

1 

i 


1 

1 

I 

I, 

li 

i: 

li 

li 

I! 

D 


where  nsnumber  of  data  points 
X.  = the  data  point. 

All  available  RMS  noise  values  for  each  day  were  averaged  to  provide  a more 
stable  RMS  noise  estimate  and  converted  to  millimicrons.  The  resulting  RMS 
noise  values  are  plotted  against  Julian  day  in  Figures  IV-1  to  lV-3  and  the 
monthly  RMS  noise  trends  derived  from  these  values  are  plotted  in  Figures 
IV-4  to  IV- 6. 

The  ANMO  short-period  RMS  noise  shows  little  change  with 
time.  The  higher  values  between  days  240  and  250  are  probably  due  to  the 
1.7  Hz  waveform  described  in  Section  III.  When  the  noise  samples  were  sel- 
ected, those  which  obviously  contained  this  waveform  were  rejected.  Appar- 
ently, some  accepted  samples  contained  a significant  amount  of  this  wave- 
form. 

The  GUMO  short-period  noise  values  up  to  approximately  day 
115  reflect  the  almost  continuous  data  clipping  caused  by  the  2000  computer 
count  per  millimicron  quantization.  Starting  at  approximately  day  125,  the 
data  was  quantized  at  2 computer  counts  per  millimicron  and  valid  estimates 
of  GUMO  RMS  noise  appear.  These  values  show  that  the  noise  field  changes 
rapidly  from  day  to  day,  maintaining  high  levels  throughout.  The  MAIO 
shor.'i-period  RMS  noise  values,  like  the  ANMO  values,  show  little  trend 
with  time.  However,  they  do  show  greater  day-to-day  variation.  The  NWAO 
and  SNZO  RMS  noise  values  both  show  an  upward  trend  with  increasing  time, 
while  the  TATO  noise  values  cover  too  limited  a time  frame  to  reveal  any  clear 
trend. 

Table  IV-1  lists  the  mean  short-period  RMS  noise  values. 

These  data  reveal  a correlation  between  the  mean  RMS  noise  level  at  a station 
and  the  distance  between  that  station  and  the  nearest  coastline.  The  lowest 
RMS  noise  values  are  recorded  at  ANMO  and  MAIO,  which  are  the  most 
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GUMO  SHORT -PERIOD  RMS  NOISE 


NWAO  SHORT-PERIOD  RMS  NOISE 
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TATO  AND  SNZO  SHORT -PERIOD  RMS  NOISE  TRENDS 


TABLE  IV- 1 

MEAN  SHORT-PERIOD  RMS  NOISE  (VERTICAL  COMPONENT) 

ALL  VALUES  IN  MILLIMICRONS  MEASURED  IN  0.  5-4.  0 Hz  PASSBAND 
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distant  of  the  evaluated  stations  from  the  nearest  coastline.  Next  in  order 
of  increasing  RMS  noise  and  decreasing  distance  to  the  nearest  coastline  is 
NWAO,  followed  by  TATO.  Although  both  GUMO  and  SNZO  are  essentially 
at  the  coastline,  their  mean  short-period  RMS  noise  values  are  substantially 
different.  (Both  are  higher  than  the  other  stations. ) This  difference  at  these 
two  stations  may  be  due  to  the  fact  that  the  coastline  nearest  SNZO  is  that  of 
the  Cook  inlet,  which  is  somewhat  sheltered  from  the  open  ocean,  while  GUMO 
lacks  any  such  shelter  from  the  oceaui. 

This  correlation  between  meam  short- period  RMS  noise  and  sta- 
tion-coastline separation  suggests  that,  above  an  RMS  'baseline'  of  about  0.  5 
millimicrons,  short-period  noise  is  predominantly  due  to  ocean  wave  energy 
injected  to  the  land  at  the  coastline. 

During  the  course  of  the  evaluation,  it  was  requested  that 
peak  noise  amplitudes  be  measured  at  one- second  period  for  the  short-per- 
iod data.  To  perform  this  task,  it  was  necessary  to  add  a module  to  the  para- 
meter measurement  section  of  the  Interactive  Seismic  Processing  System  on 
the  PDP-15  used  in  short-period  data  analysis.  Once  this  was  programmed, 
installed,  and  tested,  measurement  of  short-period  noise  amplitudes  was 
initiated.  The  noise  amplitudes  were  measured  as  maximum  zero-to-peak 
displacement  of  1- second  period  waveforms  in  the  gate  preceding  the  signal 
waveform  detected  by  the  short-period  automatic  detector.  This  noise  gate 
was  normally  approximately  20  seconds  in  length.  The  statistics  of  these 
measurements  are  reported  in  Table  IV-2  in  terms  of  the  mean  and  standard 
deviation  of  the  logarithm  of  the  measured  values.  Note  that  exponentiating 
these  means  yields  the  geometric  mean  of  the  data. 

C.  THREE-COMPONENT  LONG- PERIOD  NOISE 

The  goals  of  the  long-period  noise  analysis  were  to  estimate 
the  long-period  RMS  noise  levels,  the  peak  25  second  noise  amplitudes,  and 
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TABLE  IV-2 
SHORT -PERIOD  NOISE 

LOGjq  [peak  one-second  noise  AMPUTUDEl  STATISTICS 


Station 

Mean 

LogjQ  (AMP) 

Standard  Deviation 
Logj^j  (AMP) 

Ntmiber  of  Samples 

ANMO 

-0.  05740 

0.16278 

51 

GUMO 

1.85473 

0. 14499 

27 

MAIO 

0. 05679 

0.15088 

75 

NWAO 

0. 95440 

0.19597 

38 

TATO 

1.49354 

0.16769 

27 

SNZO 

1.79363 

0. 17075 

50 

1 
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the  spectral  content  of  the  noise  field  for  each  of  the  three  components  of 
motion  (V,  N,  and  E)  at  each  station  under  evaluation.  For  each  station,  the 
evaluation  time  frame  was  from  the  date  at  which  the  station  began  record- 
ing data  (or  1 January  1976,  whichever  was  later)  to  30  September  1976. 

Noise  samples  were  selected  by  reviewing  the  NORSAR  event 
bulletins  for  daily  time  intervals  during  which  no  seismic  events  would  be 
expected  to  be  recorded  at  a given  station.  A 4096- second  noise  sample 
was  then  processed  as  described  in  Section  II.  All  noise  samples  were  kept 
in  the  recorded  V,  N,  E configuration.  Plots  of  each  noise  sample  were 
examined  for  unreported  signals  and  for  system  malfunctions  such  as  spikes. 

Due  to  constraints  on  computer  and  analyst  time,  if  either  a signal  or  mal- 
function rendered  the  noise  sample  unusable,  that  noise  sample  was  discarded 
and  no  replacement  was  attempted.  If  neither  seismic  signals  nor  malfunctions 
were  observed,  the  entire  4096- second  sample  was  included  in  the  noise  data 
base.  In  some  cases,  only  part  of  the  noise  sample  was  unusable  due  to  a sig- 
nal or  malfunction.  When  this  occurred,  2048  seconds  of  the  sample  were 
included  in  the  noise  data  base. 

Those  noise  samples  which  passed  this  visual  inspection  were 
then  input  to  a program  which  performs  the  following  functions: 

• Compute  RMS  noise  values  uncorrected  for  instrument  response 
in  the  0.  023-0.  059  Hz  passband 

• Compute  RMS  noise  values  corrected  for  instrument  response 
using  Parseval's  formula  in  the  0.  023-0.  059  Hz  passband 

• Measure  zero-to-peak  25- second  noise  amplitudes 

• Compute  the  power  spectrum  for  each  trace  of  the  noise 
sample  and  smooth  to  128  frequencies. 
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All  measured  and  computed  values  are  output  on  punched  cards 
to  facilitate  input  to  succeeding  plot  programs.  The  25- second  noise  ampli- 
tudes are  measured  as  the  maximum  25- second  noise  amplitude  of  the  entire 
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noise  sample. 

The  long-period  RMS  noise  values  in  millimicrons  (uncorrected 
for  instrument  response)  for  the  vertical,  north,  and  east  components  of  mo- 
tion are  plotted  versus  Julian  day  in  Figures  IV-7  to  IV-12.  These  values 
were  measured  in  the  0.  023-0.  059  Hz  passband.  In  these  figures,  values  for 
consecutive  days  are  connected  by  solid  lines.  Gaps  of  one  or  two  days  in 
length  are  generally  due  to  bad  noise  samples.  Longer  gaps  are  generally 
due  to  down-time  at  the  station,  resulting  in  data  not  being  available  for 
analysis.  The  mean  values  and  associated  standard  deviations  for  these 
RMS  noise  values  are  listed  in  Table  IV- 3. 


Figures  IV-7,  IV-8,  and  IV-9  for  AMMO,  GUMO,  and  MAIO 
cover  the  samb  time  period.  For  ANMO  and  MAIO,  the  great  majority  of 
the  larger  excursions  from  the  mean  occur  during  the  days  1 to  90.  In  con- 
trast, at  GUMO  the  larger  excursions  from  the  mean  are  scattered  through- 
out the  time  period  covered.  The  implication  of  this  observation  is  that  cold 
weather  and  winter  storms  increase  the  noise  level  at  ANMO  and  MAIO.  Since 
GUMO  is  located  much  farther  south  than  these  two  stations,  it  does  not  suffer 
from  this  problem. 

Since  tiie  stations  were  not  operational,  no  noise  data  exists 
for  stations  NWAO,  TATO,  and  SNZO  for  the  early  months  of  1976.  Figures 
IV- 10,  rV-11,  and  IV-12  for  these  stations  shows  no  particular  pattern  in  the 
RMS  noise  level  excursions  from  the  mean  for  the  montiis  in  which  data  were 
available. 

Table  IV-3  presents  the  mean  long-period  RMS  noise  amplitude 
lor  the  stations  under  evaluation.  The  long- period  noise  does  not  *show  the 
strong  correlation  between  noise  level  and  station- nearest  coastline  separation 
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TABLE  IV- 3 

MEAN  LONG- PERIOD  RMS  NOISE  AMPLITUDES  IN  MILLIMICRCWS 


Station 

Mean  RMS  Noise  Amplitudes 

No.  of 
Samples 

Vertical 

North 

East 

Mean 

S.  D. 

Mean 

S.  D. 

Mean 

S.  D. 

ANMO 

9.  56 

3.83 

8.  66 

4.65 

9.  76 

3.87 

221 

GUMO 

11. 25 

3.45 

17.  20 

5.34 

18.23 

5.  34 

151 

MAIO 

8.20 

4.  37 

7.89 

3.45 

8.19 

4.  00 

207 

NWAO 

13.44 

4.  85 

17.  29 

5.67 

11.61 

4.79 

167 

TATO 

14.  22 

6.04 

15.65 

7.14 

18. 12 

9.  56 

95 

SNZO 

45.92 

17.  28 

27.  26 

14.  22 

30. 17 

14.41 

148 

S.  D.  - Standard  Deviation 


th*t  waa  found  for  the  short-period  noise.  The  stations  can  be  grouped  by  noise 
level,  however.  Stations  ANMO  and  MAIO  form  the  group  with  lower  noise 
levels  and  stations  GUMO,  NWAO,  and  TATO  form  the  goup  with  higher  noise 
levels.  Considering  the  locations  of  these  stations,  this  suggests  that  the 
somewhat  higher  noise  levels  at  GUMO,  NWAO,  and  TATO  are  due  to  wave 
energy  injected  at  the  coastline.  (Station  SNZO  was  not  included  in  this  group- 
ing, since  the  noise  levels  at  this  station  seem  too  high  to  be  due  to  seismic 
noise. ) 

Figures  IV- 13  to  IV- 18  show  the  monthly  noise  level  trends  for 
the  stations  under  evaluation.  Stations  ANMO  and  MAIO  show  a downward 
trend  in  noise  level  from  January  1976  to  September  1976.  It  seems  reasonable 
to  ascribe  this  trend  to  the  decrease  in  storm  activity  with  the  passing  of  the 
winter  months  and  associated  winter  storms.  Station  GUMO  shows  a general 
upward  trend  in  noise  level  during  the  same  time  period.  Considering  the 
location  of  the  station,  it  seems  possible  that  this  increase  in  seismic  noise 
is  due  to  wave  energy  generated  by  tropical  storms.  Stations  NWAO  and  SNZO 
show  a general  upward  trend  in  noise  levels  from  their  start-up  dates  to  Sept- 
ember 1976,  dropping  off  slightly  after  peaking  in  July  1976.  Since  these  sta- 
tions are  In  the  southern  hemisphere,  this  represents  an  increase  in  noise 
level  with  the  onset  of  the  winter  months.  Again,  this  increased  winter  noise 
can  be  associated  with  increased  storm  activity.  Disregarding  the  noise  level 
for  May  1976,  station  TATO  shows  an  upward  trend  in  noise  level  in  the  same 
manner  as  station  GUMO.  Thus,  this  increase  can  be  ascribed  to  increased 
tropical  storm  activity  in  the  late  summer  months. 

Figures  IV- 19  to  IV- 24  show  the  average  RMS  amplitude  spectra 
for  the  noise  at  each  station  under  evaluation.  These  spectra  are  uncorrected 
for  instrument  response.  They  were  determined  by  averaging  the  individual 
spectra  computed  from  the  noise  samples  and  converting  them  to  RMS  ampli- 
tudes using  Parseval's  formulas 
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FIGURE  IV- 13 

ANMO  LONG-PERIOD  RMS  NOISE  TRENDS 
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FIGURE  IV- 19 

AVERAGE  RMS  AMPLITUDE  SPECTRA- ANMO  LONG- PERIOD  NOISE 
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FIGURE  IV- 21 

AVERAGE  RMS  AMPLITUDE  SPECTRA-MAIO  LONG- PERIOD  NOISE 
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RMS^  = 
a 

where  if  * 

|A(f.)|^ 

C(£.) 

1 

a = 

b 

Since  the  RMS  amplitude  at  each  discrete  frequency  was  desired,  in  this  case 
a=b.  Also,  since  no  instrument  response  corrections  were  made,  C(fj)=l  for 
all  frequencies. 

The  left-hand  side  of  each  figure  shows  the  average  RMS  noise 
amplitude  spectra  for  each  component  of  motion.  These  spectra  were  computed 
for  the  period  range  of  11. 6 to  85.  3 seconds  so  that  the  analyst  can  inspect  the 
noise  levels  in  and  around  the  signal  periods.  The  right-hand  side  of  each 
figure  shows  the  log  mean  spectra  with  vertical  bars  representing  plus  and 
minus  one  standard  deviation.  These  are  presented  to  give  a qualitative 
idea  of  the  day-to-day  variability  of  the  noise  from  period  to  period. 


the  elemental  frequency  interval  ( feO.  001953  Hs) 

the  discrete  Fourier  transform  spectral  density 
estimate  at  frequency  f., 

the  instrument  response  correction  at  frequency  f., 
the  initial  frequency  index,  and 
the  final  frequency  index. 


The  following  comments  can  be  made  about  the  individual  fig- 

Figure  IV-19---ANMO  noise:  The  three  components  show 
nearly  identical  spectral  content  and  level  of  the  local  noise 
field.  The  noise  fitld  has  its  greatest  day-to-day  variability 
between  14  and  25  seconds  period. 

Figure  lV-20 GUMO  noise:  The  three  components  ahow 

markedly  different  spectral  content  and  level  of  the  local  noise 
field  from  one  component  to  the  next.  The  horisontal  components 
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show  higher  noise  levels  and  day-to-day  variability  then 
does  the  vertical  component. 

Figure  IV-21  - MAIO  noise:  At  periods  below  30  seconds,  the 
three  components  show  nearly  identical  spectral  content  and 
level  of  the  local  noise  field.  The  horizontal  components  show 
a secondary  peak  at  37  seconds  nearly  as  large  as  the  micro- 
seismic  peak.  The  noise  field  has  its  greatest  variability  be- 
tween 14  and  23  seconds  period  for  all  components  and  at  37 
seconds  period  for  the  horizontal  components. 

Figure  IV-22  - NWAO  noise:  The  three  components  of  noise 
recorded  at  this  station  show  more  variation  in  spectral  content 
and  level  of  the  noise  from  one  to  the  next  than  do  the  three  com- 
ponents of  noise  recorded  at  ANMO  or  MAIO.  The  noise  field 
has  its  greatest  day-to-day  variability  between  13  and  23  seconds 
period. 

Figure  IV-23  - TATO  noise:  At  periods  below  35  seconds,  the 
three  components  show  similar  spectral  content  and  level  of  the 
noise.  The  component-to- component  variation  appears  no  greater 
than  that  observed  at  NWAO.  Above  35  seconds,  the  horizontal 
components  show  rapidly  increasing  noise  levels  while  the  ver- 
tical component  remains  fairly  constant.  The  day-to-day  varia- 
tion appears  to  be  quite  large  at  all  frequencies. 

Figure  IV-24  - SNZO  noise:  At  periods  below  20  seconds,  the 
three  components  show  similar  spectral  content  and  level  of  the 
noise.  The  component-to-component  variation  appears  no  greater 
than  that  observed  at  NWAO  and  TATO.  At  periods  above  20  sec- 
onds, the  vertical  component  noise  increases  rapidly,  causing 
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the  high  RMS  noise  values  observed  in  Table  IV- 3.  The  day 
to-day  variation  appears  to  be  quite  large  at  all  frequencies 


The  above  comments  lead  the  author  to  believe  that  oceanic 
wave  energy  plays  a large  part  in  determining  the  level,  spectral  content,  and 
variability  of  the  local  long-period  noise  field  for  stations  near  the  coast. 
Stations  ANMO  and  GUMO  of  the  six  evaluated  show  low  noise  levels,  very 
similar  noise  spectra  for  the  three  components,  and  day-to-day  noise  varia- 
tions which  are  large  only  near  the  microseismic  peak.  As  the  station- coast- 
line separation  decreases,  the  noise  levels  rise,  the  differences  among  the 
spectra  of  the  three  components  increase,  and  the  day-to-day  variability  of 
the  noise  rises.  While  the  effect  of  oceanic  wave  energy  on  the  long-period 
noise  field  is  not  as  severe  as  it  is  on  the  short-period  noise  field,  it  is  in 
the  author's  opinion  significant  enough  to  affect  the  detection  capability  of  the 
station.  The  spectral  horizontal  components  at  longer  periods  appear  to  be 
abnormally  high  when  compared  to  those  of  other  stations.  This  is  apparently 
due  to  the  seismometer  being  underwater. 


The  level  and  spectral  content  of  the  vertical  component  of  the 
SNZO  noise  field  do  not  appear  reasonable.  The  filtering  of  the  broadband 
data  to  produce  the  long-period  data  cut  down  the  64- second  data  by  approxi- 
matelv  9 dB  relative  to  the  25- second  data.  However,  Figure  IV- 24  shows 
that  after  this  filtering,  the  64- second  data  are  approximately  7 dB  above 
the  25- second  data.  These  high  noise  levels  at  longer  periods  have  been 
ascribed  to  electronic  noise  on  the  vertical  component. 


When  the  average  noise  spectra  for  GUMO  were  plotted,  it  was 

□ 

noticed  that  they  were  quite  different  from  those  prepared  for  the  preliminary 
Seismic  Research  Observatories  evaluation,  which  were  made  by  averaging 
the  spectra  of  noise  samples  recorded  during  the  period  January  to  April, 
1976.  Figure  IV-25  was  created  to  investigate  this  difference.  Each  plot  of 
this  figure  shows  the  average  spectra  for  noise  samples  recorded  during  one 


month  of  1976,  (The  number  at  the  top  of  each  plot  indicates  the  month.  ) 
Although  the  noise  spectra  for  the  first  four  months  show  marked  variations 
from  one  to  the  next,  they  form  a group  distinctly  different  from  the  noise 
spectra  for  the  succeeding  months.  These  months  show  much  greater  varia- 
tions in  noise  level  and  spectral  content  than  do  the  first  four  months.  Thus, 
it  is  not  surprising  that  the  nine-month  average  presented  in  this  report  dif- 
fers so  much  from  the  four-month  average  presented  last  year. 

Peak  25-second  noise  amplitudes  were  measured  on  all  noise 
samples  during  the  course  of  the  noise  analysis.  Each  value  measured  rep- 
resents the  absolute  value  of  the  largest  zero-to-peak  excursion  of  a wave- 
form of  period  25  seconds  plus-or-minus  one  second.  Table  IV-4  lists  the 
mean  and  standard  deviation  of  the  logarithms  of  the  measured  peak  ampli- 
tudes for  each  station.  They  are  presented  in  this  form  since  it  is  in  this 
form  that  they  are  used  in  the  evaluation  of  the  Seismic  Research  Observatory 
network  (Section  VII). 
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SECTION  V 

SRO  DETECTION  CAPABILITY 
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A.  DISCUSSION 

In  paat  evaluation  task*,  detection  statistics  were  derived 
from  an  event  population  for  every  member  of  which  a clear  detection/ 
non-detection  decision  could  be  made;  i,  e. , the  analyst  could  state  that  he 
either  saw  the  sought-for  signal  or  he  saw  seismic  noise.  Unfortunately, 
the  world  does  not  always  present  the  analyst  with  such  a clear-cut  case. 
Mixed  events,  system  failures  resulting  in  no  data  being  recorded,  and 
malfunctions  all  tend  to  obscure  the  detection  capability  picture. 

The  term  'mixed  event*  refers  to  the  case  where  the  sought- 
for  signal  is  obscured  or  completely  masked  by  a second  signal.  This  can 
happen  either  when  the  two  signals  arrive  at  the  station  at  essentially  the 
same  time  or  when  a larger  signal  arrives  before  the  signal  under  analysis, 
burying  this  signal  in  its  coda.  The  term  'system  failure'  refers  to  the 
total  shutting  down  of  the  station  so  that  no  data  is  recorded.  The  term 
'malfunction'  refers  to  any  partial  failure  of  the  system,  from  sensor  unit 
to  reception  of  data  at  the  Seismic  Data  Analysis  Center,  which  causes  de- 
gradation of  the  seismic  data. 

The  problem  of  mixed  events  is  often  difficult  for  the  analyst 
to  resolve  and  is  probably  the  nutjor  source  of  false  alarms.  (The  term 
'false  alarm'  in  this  case  means  declaring  a detection  when  in  fact  the  ob- 
served signal  is  from  an  event  other  than  that  under  analysis. ) When  a sig- 
nal is  observed  in  the  time  gate  of  the  event  under  analysis,  the  analyst  first 
checks  the  waveforms  on  the  three  components  of  motion  to  see  that  their 


V-1 


inter-relationship*  are  correct.  If  still  in  doubt,  the  analyst  checks  avail- 
able event  lists  to  see  whether  any  other  reported  event  could  have  arrived 
in  the  signal  gate.  In  general,  the  analyst  declares  a detection  if  a dispersed 
signal  is  observed  having  the  correct  inter-relationships  between  the  L>ove 
and  Rayleigh  waves  and  if  no  other  event  has  been  reported  which  could  be 
mistaken  for  the  event  under  analysis. 

Since  mixed  events  and  (secondarily)  events  for  which  no 
data  were  recorded  or  which  contained  malfunctions  are  a fairly  common 
problem,  the  Seismic  Research  Observatory  detection  capability  estimates 
are  calculated  in  two  ways.  The  first  of  these  is  termed  the  absolute  de- 
tection capability  estimate.  When  computing  this  estimate,  all  mixed  events, 
events  for  which  no  data  were  recorded,  and  events  containing  malfunctions 
are  counted  as  non-detections  when  forming  the  detection  statistics.  This 
approach  gives  a real-world  detection  capability  estimate. 

The  second  of  these  estimates  is  termed  the  conditional  de- 
tection capability  estimate.  When  computing  this  estimate,  all  mixed  events, 
events  for  which  no  data  were  recorded,  and  events  containing  malfunctions 
are  rejected  from  the  detection  statistics.  This  approach  gives  an  ideal  de- 
tection capability  estimate.  The  value  of  this  ideal  estimate  is  that  it  shows 
the  detection  capability  improvement  possible  if  the  reliability  of  the  instru- 
mentation can  be  improved  and  if  methods  of  separating  mixed  events  can 
be  found. 

The  number  used  to  represent  detection  capability  is  the 
50  percent  detection  threshold,  denoted  by  The  50  percent  detection 

threshold  is  the  body  wave  magnitude  for  which  the  probability  of  detection  is 
0.  5.  It  is  computed  by  fitting  the  Gaussian  probability  function  to  the  detection 
statistics  by  maximum  likelihood  method  (Ringdal,  1974). 
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It  must  be  kept  in  mind  that,  since  the  data  base  was  derived 
from  the  NORSAR  event  bulletins,  all  detection  capability  estimates  are  in 
terms  of  NORSAR  m^  units.  This  fact  is  impKjrtant,  since  NORSAR  m^  units 
are  not  the  same  as  those  from  other  earthquake  parameters  sources  such  as 
the  National  Earthquake  Information  Service  (NEIS)  event  list. 

Depending  on  the  quantity  of  data  available,  one  of  three  levels 
of  confidence  is  placed  on  each  50  percent  detection  threshold.  If  a 50  percent 
detection  threshold  is  given  to  two  decimal  places  (e.  g. , m = 4.  56), 
sufficient  detection  statistics  were  available  to  compute  a reliable  detection 
capability  estimate.  If  a 50  percent  detection  threshold  is  enclosed  in  paren- 
theses and  given  to  only  one  decimal  place,  (e.  g.  , ^))  detection 

statistics  were  sparse  and  the  detection  capability  estimate  can  be  considered 
to  be  only  a first  approximation.  Finally,  if  a 50  percent  detection  threshold 
is  expressed  as  greater  than  some  value  (e.  g. , m^gQ>5.  0),  the  detection 
statistics  were  too  sparse  to  allow  fitting  of  the  Gaussian  probability  curve. 

In  most  cases,  this  is  due  to  a lack  of  detected  events.  When  this  occurred, 
the  value  given  is  the  m^^  of  the  largest  non-detected  event. 

t 

B.  SHORT-PERIOD  DETECTION  CAPABILITY  ESTIMATES  j 

The  criteria  for  determining  whether  an  event  was  detected  on  I 

the  recorded  vertical  short-period  component  are  as  follows:  J 

IT 

• The  waveform  is  at  least  12  dB  above  the  surrounding  noise 

waveforms  ^ 

• The  waveform  is  impulsive,  not  emergent  ^ 


• The  waveform  begins  within  ± 20  seconds  of  the  predicted 

arrival  time. 

The  requirement  that  the  maximum  waveform  amplitude  be 
at  least  12  dB  above  the  surrounding  noise  waveforms  is  set  high  because 
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the  amplitude  change  is  of  overriding  importance  when  searching  for  P- 
waves  on  a single  component  of  motion.  It  must  be  remembered  that  the 
analyst  does  not  have  available  the  information  which  three  components  of 
short-period  motion  provide,  i.  e.  , the  approximate  arrival  azimuth  pro- 
vided by  the  direction  of  motion  of  the  first  arrival  and  the  approximate  epi- 
central  distance  provided  by  the  P-S  arrival  time  separation.  (S  is  rarely 
observed  on  the  vertical  component. ) Therefore,  the  analyst  must  rely  on 
an  abrupt  amplitude  change  to  feel  confident  that  the  waveform  under  consid- 
eration is  a P-wave.  Occasionally,  it  is  possible  to  pick  a P-wave  when  this 
criterion  is  not  satisfied.  This  occurs  when  a change  in  frequency  and  am- 
plitude are  simultaneously  observed.  Typically,  this  change  in  frequency  is 
from  two  Hertz  to  one  Hertz. 

The  impulsive  waveform  requirement  is  made  in  order  that  the 
analyst  be  able  to  determine  the  P-wave  arrival  time.  Without  this,  the  an- 
alyst cannot  be  confident  of  fulfilling  the  third  criterion.  This  third  criterion, 
the  £ 20  second  P-wave  arrival  gate,  is  an  implicit  recognition  that  the  origin 
time  and/or  epicenter  coordinates  may  be  in  error,  resulting  in  a deviation  of 
the  actual  arrival  time  from  the  predicted  arrival  time.  Also,  it  must  be  noted 
that  the  predicted  arrival  time  is  based  on  the  assumption  that  the  event  occur- 
red at  normal  (33  km)  depth.  Departures  from  this  normal  depth  will  also 
cause  departures  from  the  predicted  arrival  time. 

The  third  detection  criterion  is  not  absolute.  It  is  possible  that 
errors  in  computation  of  origin  time,  location,  and  depth  may  all  combine  to 
pla  ce  the  actual  P-wave  arrival  outside  the  £ 20  second  gate.  In  cases  where 
a P-wave  was  observed  outside  the  gate  and  no  other  event  could  be  found  in 
the  NORSAR  bulletins  which  would  have  arrived  at  that  time,  a detection  was 
declared.  This  was,  however,  a rare  occurrence. 

Short-period  50  percent  detection  thresholds  were  first  com- 
puted by  station  for  all  regions  as  a unit  using  tiuree  sub- sets  of  the  detection 
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statistics.  The  first  sub-set  ('ALL  REGIONS  A')  consists  of  all  short- period 
events  with  epicentral  distances  less  than  103  degrees.  This  produces  the 
previously  defined  absolute  detection  capability  estimate  for  each  station.  The 
second  sub- set  ('ALL  REGIONS  B')  produces  the  previously  defined  conditional 
detection  capability  estimate.  Note  that  this  estimate  contains  the  non-detec- 
tions caused  by  the  failure  of  the  automatic  detector  to  trigger.  These  de- 
tection capability  estimates  represent  the  short-period  detection  capability 
of  the  stations  as  they  are  currently  constituted.  The  third  sub- set  ('ALL 
REGIONS  C')  contains  only  those  short- period  events  for  which  data  were 
recorded.  (The  non-detections  in  this  subset  were  found  primarily  in  'false 
alarm'  records  output  by  the  automatic  detector.  ) Tnese  detection  capability 
estimates  were  calculated  in  an  attempt  to  estimate  the  effect  of  the  automatic 
detector  on  detection  capability.  The  detection  capability  estimates  for  the 
individual  regions  were  computed  in  the  same  fashion  as  they  were  for  'ALL 
REGIONS  B',  i.e. , those  events  not  detected  by  the  automatic  detector  were 
counted  as  non-detections  in  the  detection  statistics. 

The  detection  statistics  and  fitted  maximum  likelihood  curves 
for  the  short-period  data  are  shown  in  Figures  V-1  to  V-8.  Each  of  the  'sub- 
figures' in  Figures  V-1  to  V-8  represents  one  short-period  detection  capability 
estimate.  The  upper  portion  of  each  'sub-figure'  is  a histogram  giving  the 
detection  statistics  as  a function  of  body-wave  magnitude  for  that  particular 
case.  The  lower  portion  shows  the  percentage  of  events  detected  at  each 
body-wave  magnitude  (represented  by  asterisks),  the  fitted  maximum  like- 
lihood curve  (represented  by  a solid  line),  and  the  90  percent  confidence  limits 
for  this  curve  (i;epresented  by  dashed  lines).  The  values  for  'MB50'  and 
'MB90'  shown  on  the  figures  are  the  50  and  90  percent  detection  thresholds  as 
picked  from  the  maximum  likelihood  curve.  The  value  shown  for  'SIGMA' 
is  the  standard  deviation  of  Gaussian  probability  function  obtained  by  maximum 
likelihood  method. 
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FIGURE  V-1 

ANMO  SHORT-PERIOD  DETECTION  STATISTICS 
ALL  REGIONS  A.  B,  AND  C 
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The  titles  for  each  sub -figure  define  the  station  under  evaluation 
and  the  region  of  the  events  from  which  the  detection  statistics  were  derived. 
The  short-period  detection  capability  estimates  derived  from  Figures  V-1  to 
V-8  are  summarized  in  Table  V-1.  In  this  table,  '^^103'  indicates  that  the 
mean  epicentral  distance  from  region  to  station  was  greater  than  103  degrees 
and  that  therefore  no  P-wave  could  be  detected.  As  noted  in  Table  11-4,  only 
two  short-period  events  were  detected  at  GUMO  and  none  at  SNZO.  This 
paucity  of  detections  can  be  attributed  to  the  high  noise  levels  at  these  stations. 

At  ANMO,  the  absolute  and  conditional  detection  capability 
estimates  differ  by  0.  2 units.  One  can  see  from  Table  n-4  that  this 
difference  is  almost  entirely  due  to  events  counted  as  non-detections  due  to  no 
data  being  recorded.  Thus,  improving  the  station  reliability  will  bring  the 
absolute  and  conditional  detection  capabilities  very  close  together.  In 
comparison,  at  MAIO  the  absolute  and  conditional  detection  capability  estim- 
ates differ  by  only  0. 1 m . units.  This  difference  is  primarily  due  to  malfunc- 

D 

tions. 

The  effect  of  the  automatic  detector  can  be  estimated  by  the 
difference  between  'ALL  REGIONS  B*  and  'ALL  REGIONS  C'.  At  ANMO  this 
difference  is  approximately  0.  5 units,  while  at  MAIO  it  is  approximately 
0.  9 units.  Part  of  this  difference  may  be  due  to  higher  noise  levels  mask- 
ing some  short-period  arrivals  (i.  e. , higher  noise  levels  vould  prevent  the 
automatic  detector  from  'seeing'  an  arrival).  However,  it  is  believed  that 
part  of  this  difference  is  due  to  improper  functioning  of  the  automatic  detector 
itself  caused  by  setting  the  tum-on  threshold  too  high.  This  can  only  be  re- 
solved by  obtaining  continuous  short-period  data,  visually  examining  this  data 
for  detections,  and  finally  running  this  data  through  the  automatic  detector  to 
determine  which  events  it  can  detect. 

The  individual  region  detection  capability  estimates  of  Table 
V-1  are  determined  in  the  same  fashion  as  the  estimates  for  'ALL  REGIONS  B'. 
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TABLE  V-1 

SRO  SP  DETECTION  CAPABILITY 


I 


Region 


111  Region*  A 
l11  Region*  B 
ill  Region*  C 
1 
2 


5.  07 
4.85 
4.  37 
4.  90 
4.  60 


>6.0 

>6.0 

>6.0 

>6.4 

>5.6 


50  Percent  Detection  Threshold 


MAIO  NWAO 


3.74  (5.1) 

5.22  >5.7 

5.04  >5.6 


5.7 
5.  5 
(4.  9)  >5.  7 

>5.7  >5.7 

>5.  0 >5.  6 


4>103  >5.7 

.A>103  >6.0 


^>103 

A>103 

(5.7) 

>4.6 


>6.0 

>6.5 


»4. 6 ^03 

(4.8)  I 


(4.  0) 
4. 16 

(4.2) 
(4.0) 

(4.3) 


>6.2 
>5.3 
>5.2 
>6.0 
.A  >103 


(4.9)  ^>103 

(4.4)  (5.4) 


(5.4)  >6.2 

(5.3)  4>103 


>5.3 

(5.5) 

>5.3 


1>103 

1>103 

1>103 


4>103 

(5.3)  >6.5 
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All  Region*  A Mixed,  no  data  recorded,  no  short-period  data 
recorded,  and  malfunction*  are  counted  a*  non- 
detection*. 

All  Region*  B Mixed,  no  data  recorded,  and  malfunction*  re- 
jected from  detection  statistic* 

All  Region*  C Mixed,  no  data  recorded,  no  short-period  data 

recorded,  and  malfunction*  rejected  from  detection 
statistic*. 

Individual  region  estimate*  were  computed  in  tiie  same  fashion  a* 

'All  Region*  B". 

50  percent  detection  threshold*  are  in  NORSAR  mj^  unit*. 
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Thu«,  these  detection  capability  estimates  represent  the  detection  capability  of 
the  stations  as  they  now  fvinction.  With  few  exceptions,  the  detection  statis- 
tics were  too  sparse  to  permit  determination  of  reliable  detection  capability 
estimates.  If  good  estimates  are  desired,  a great  deal  more  short-period 
data  (at  most  stations,  double  or  more  the  amount  already  processed)  must  be 
processed  and  analysed. 

C.  LONG-PERIOD  DETECTION  CAPABILITY  ESTIMATES 

Estimates  of  long-period  Seismic  Research  Observatory  de- 
tection capability  are  presented  for  each  of  the  six  stations  evaluated  in  this 
report.  The  criteria  for  determining  whether  detection  has  been  achieved 
for  a given  event  are: 

• The  presence  of  dispersion  in  the  signal  gate 

a The  presence  of  a peak  in  the  dispersed  wave  train  3 dB  or 

more  above  any  peak  outside  the  dispersed  wave  train  and 
inside  a time  gate  starting  600  seconds  before  the  predicted 
Love  wave  arrival  time  and  ending  600  seconds  after  the 
estimated  Rayleigh  wave  end  time. 

• The  occurrence  of  signal  onset  within  ^180  seconds  of  the 
predicted  signal  onset  time 

• Detection  of  the  event  on  at  least  two  components. 

Occasionally,  an  event  is  considered  to  be  detected  which  does  not  satisfy 
all  three.  Signal  peaks  are  occasionally  less  than  3 dB  above  the  noise  peaks 
and  still  recognized  as  signals  from  their  dispersion  characteristics.  Also, 
at  a given  station  it  is  sometimes  possible  to  find  specific  features  of  a 
seismic  watreform  from  a given  region,  which  enable  the  analyst  to  detect 
the  event  even  though  not  all  detection  criteria  have  been  satisfied.  An 
interesting  example  of  this  is  discussed  in  the  preliminary  Seismic  Research 
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Observatory  report  (Strauss,  1976).  The  last  criterion  is  intended  to  reduce 
the  number  of  false  alarms,  i.  e. , calling  an  event  detected  when  it  is  not. 


I The  long-period  detection  statistics  and  derived  maximum 

likelihood  curves  are  presented  in  Figures  V-9  to  V-24.  The  50  percent 
I detection  thresholds  were  first  computed  by  station  for  all  regions  as  a unit 

using  two  subsets  of  the  detection  statistics.  The  first  subset  ('ALL  REGIONS 
I A')  is  used  to  estimate  the  previously  defined  absolute  detection  capability 

estimate.  In  this  case,  mixed  events,  events  for  which  no  data  were  recorded, 
I and  events  for  which  the  data  contained  malfunctions  were  counted  as  non- 

^ detections.  The  second  subset  ('ALL  REGIONS  B')  is  used  to  estimate  the  pre- 

I viously  defined  conditional  detection  capability  estimate.  In  this  case,  mixed 

events,  events  for  which  no  data  were  recorded,  and  events  for  which  the  data 
I contained  malfunctions  were  rejected  from  the  detection  statistics.  The 

estimates  for  the  individual  regions  are  conditional  detection  capability  es- 
i timates  deternnined  in  the  same  fashion  as  'ALL  REGIONS  B'. 
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Figures  V-9  to  V-24  are  laid  out  in  the  same  manner  as  Fig- 
ures V-1  to  V-8  for  the  short-period  detection  capability  estimates.  Each 
of  the  'sub-figures'  represents  one  long-period  detection  capability  estimate, 
where  the  upper  portion  shows  a histogram  of  the  detection  statistics  and  the 
lower  portion  shows  the  maximum  likelihood  curve  fitted  to  these  statistics. 

The  long-period  detection  capability  estimates  derived  from  these 
figures  are  summarized  in  Table  V-2  . 


The  differences  between  the  absolute  and  conditional  detection 
capability  estimates  fall  into  three  groups.  At  NWAO  and  SNZO,  this  differ- 
ence is  approximately  0.  3 m^  units.  At  NWAO,  this  is  primarily  due  to 
mixed  events,  while  at  SNZO  system  failures  resulting  in  no  data  being  re- 
corded and  malfunctions  have  roughly  as  much  effect  as  do  ^e  mixed  events. 
(See  Table  II-4.)  At  ANMO,  MAIO,  and  TATO,  the  difference  between  the 
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FIGURE  V-10 

ANMO  LONG- PERIOD  DETECTION  STATISTICS 
REGIONS  2,  5,  6,  AND  7 


V-19 


FIGURE  V-16 


MAIO  LONG- PERIOD  DETECTION  STATISTICS 
REGIONS  2,  5,  6,  AND  7 


FIGURE  V-17 


MAIO  LONG- PERIOD  DETECTION  STATISTICS 
REGIONS  9.  10,  11.  AND  12 


TABLE  V-2 

SRO  LONG-PERIOD  DETECTION  CAPABILITY 


Region 

ANMO 

GUMO 

MAIO 

NWAO 

TATO 

SNZO 

All  Regions  A 

S.  04 

5.09 

4.55 

5.20 

4.  50 

All  Regions  B 

4.  61 

4.53 

4.  13 

4.  91 

4.12 

1 

4. 82 

4.  57 

4.48 

(5.5) 

(4.3) 

(5.2) 

2 

(4.  8) 

4.24 

4.39 

(5.0) 

(4.1) 

(4.8) 

5 

(4.1) 

(4.2) 

(4.2) 

(4.4) 

(3.9) 

>4.4 

6 

(4.  5) 

(4.5) 

3.80 

4.  57 

(3.9) 

>4.  9 

7 

(4.7) 

(4.5) 

(3.8) 

(5.1) 

(4.1) 

>5.2 

9 

4.47 

(5.0) 

3.61 

(5.2) 

(4.8) 

>5.2 

10 

4.25 

(4.  6) 

3.  98 

(4.7) 

(4.3) 

>4.7 

11 

(4.2) 

>4.3 

3.97 

>5.0 

4.0 

>4.5 

12 

(4.3) 

>4.3 

(3.8) 

>4.5 

(3.6) 

>4.9 

All  Regiona  A Mixed  events,  eventa  for  which  no  data  were  re-  * t 

corded  and  events  containing  malfunctions  are  | 

counted  as  non-detections. 

All  Regions  B Mixed  events,  events  for  which  no  data  were  re-  | 

corded,  and  events  containing  malfunctions  are  ^ > 

rejected  from  the  detection  statistic.  ^ 

50  percent  detection  thresholds  are  in  NORSAR  m^  units.  - ^ 

lI 


absolute  and  conditional  detection  capabilities  is  approximately  0.4  m^  units. 
In  each  case,  this  difference  is  due  more  to  mixed  events  than  to  system 
failures  or  malfunctions.  At  GUMO,  the  difference  between  the  absolute  and 
conditional  detection  capability  estimates  is  greater  than  0.  5 units.  The 
increase  in  this  difference  over  that  of  the  previous  group  of  stations  is  due 
to  the  increased  system  failure  rate  (resulting  in  no  data  being  recorded) 
of  GUMO.  The  remainder  of  the  difference  is  primarily  due  to  mixed  events. 
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EARTHQUAKE- PRESUMED  EXPLOSION  DISCRIMINATION 
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A.  DISCUSSION 

I 

This  section  considers  the  question  of  discriminating  between 
earthquakes  and  presumed  nuclear  explosions  using  long-period  data.  The  dis- 
crimination method  used  is  the  surface  wave  magnitude  (M^)  versus  body  wave 
magnitude  (m^)  plot.  This  plot  of  M^  versus  m^  is  expected  to  function  as  an 
earthquake -presumed  nuclear  explosion  discriminant  since,  for  a given  m^, 
an  explosion  produces  much  lower  Rayleigh  and  Love  waves  than  does  an  earth- 
quake. (Theoretically,  an  explosion  should  generate  no  Love  waves,  since 
the  source  is  completely  compressional.  However,  some  Love  wave  energy 
is  radiated  from  tectonic  strain  release  (Sun,  1977)  subsequent  to  the  explosion.  ) 
Therefore,  a plot  of  M versus  m^  can  be  expected  to  show  a separation  of  the 
data  points  into  an  earthquake  population  and  an  explosion  population. 

The  data  base  for  this  discrimination  study  is  comprised  of 
all  events  which  were  visually  detected  on  bandpass-filtered  (0.023-0.  059  Hz 
passband)  plots. 

B.  COMPUTATION  OF  SURFACE  WAVE  MAGNITUDES 

In  earlier  evaluation  tasks,  the  signal  amplitude  and  period 
values  used  in  computing  surface  wave  magnitudes  were  measured  man- 

ually on  the  filtered  signal  plots  at  various  periods.  The  computation  of  M^ 
then  used  the  equation: 
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where: 

A = peak-to-peak  amplitude  measured  in  inches  on  the 
plot 

SF  = plot  scale  factor  in  computer  counts  per  inch 
T = period  in  seconds  of  the  measured  amplitude 
Q = quantization  factor  = 5 computer  counts  per  millimicron 

G = instrument  response  correction  factor,  and 

L a epicentral  distance  in  degrees. 

This  approach  has  several  disadvantages.  First,  the  analyst  is 
apt  to  have  difficulties  finding  the  largest  peak  at  each  desired  period- -it  is 
sometimes  necessary  to  measure  the  periods  of  a number  of  waveforms  before 
finding  the  desired  period.  Next,  having  found  the  waveform  for  which 
is  to  be  measured,  the  analyst  may  make  a measurement  error.  Finally, 
errors  are  apt  to  occur  either  in  transcribing  the  measured  values  or  in  cal- 
culating the  surface-wave  magnitude  from  these  values. 

In  order  to  avoid  these  problems  and  make  the  process  of  ob- 
M,  vUu..  1...  . prog,™  w..  written  to  perform  the  me..- 

urements  and  calculations  automatically.  This  program  operates  by  first 
finding  the  times  of  all  zero  crossings  in  a given  time  gate  and  the  maximum 
absolute  amplitude  in  millimicrons  between  each  pair  of  adjacent  zero  cross- 
ings. The  waveform  period  is  then  simply 


T = 2 ♦ 


Time  of 


Zero. - Time 
i+l 


of  Zero^^ 


The  program  reads  the  epicentral  distance  from  the  event  header.  The  quan- 
tization factor  and  instrument  response  corrections  are  built  into  the  program. 
(The  instrument  response  corrections  are  derived  from  the  instrument  response 
curve  of  Figure  VI-1. ) 

The  program  uses  the  above  information  to  compute  a surface- 
wave  magnitude  for  each  amplitude -period  pair.  These  values  with  their 
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NORMALIZED  INSTRUMENT  RESPONSE  CURVES 
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associated  periods  and  times  of  measurement  are  printed  out.  Finally  the 
program  picks  the  largest  value  at  each  specified  period  of  interest  and 
prints  it  out.  The  analyst  now  checks  these  largest  values  against  the 
plot  of  the  data.  If  he  does  not  like  an  value  at  a particular  period  (for 

example,  the  M value  appears  to  be  associated  idth  a noise  pulse  in  the 

s 

signal  gate)  he  can  select  another  value  from  the  list  generated  by  the 
program. 

By  using  this  program,  measurement  errors  have  been  elim- 
inated and  transcription  errors  greatly  reduced  (since  the  analyst  only  writes 
down  the  final  value  and  none  of  the  intermediate  computational  values). 
Also,  the  time  required  to  arrive  at  the  values  for  an  event  has  been  re- 
duced. 

A further  benefit  to  the  evaluation  task  was  found  once  this 
program  was  put  into  use.  Since  measurements  were  no  longer  being  made 
directly  on  the  plot,  it  was  possible  to  reduce  the  length  of  the  plots  by  one- 
half.  When  measurements  were  made  on  the  plot,  a horizontal  scale  of  100 
seconds  per  inch  was  considered  necessary  to  minimize  measurement  errors. 
Plots  are  now  made  with  a horizontal  scale  of  200  seconds  per  inch.  This 
reduced  the  amount  of  plot  paper  and  plot  time  by  one-half  and  produced  more 
manageable  plots. 

During  the  preliminary  Seismic  Research  Observatories  eval- 
uation, was  measured  at  periods  of  20,  30,  and  40  seconds.  However, 
measurable  40  second  energy  was  not  often  found- -only  half  as  many  40  second 

M values  were  measured  as  20  second  or  30  second  M values.  It  was  also 

s » 

noted  that  events  close  to  MAIO  have  the  majority  of  their  energy  concentrated 
at  or  near  25  seconds  period.  For  these  reasons,  it  has  been  decided  to  drop 
the  measurement  of  40  second  M and  add  the  measurement  of  25  second  M . 
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DISCRIMINATION  RESULTS 


Table  VI- 1 lists  those  events  of  the  data  base  which  may  be 
termed  presumed  nuclear  explosions.  They  were  selected  on  the  basis  of 
their  epicentral  locations  and  body  wave  magnitudes.  Thus,  the  events  from 
Region  16  were  selected  because  their  epicenters  are  at  the  Nevada  Test 
Site  and  their  body  wave  magnitudes  are  much  larger  than  known  earthquakes 
in  the  area.  The  events  from  Region  8 were  selected  because  their  epicenters 
are  at  the  eastern  Kazakh  test  site  and  their  body  wave  magnitudes  are  larger 
than  most  earthquakes  from  that  area.  The  one  event  from  Region  10  has  its 
epicenter  in  an  area  of  low  seismicity  where  peaceful  nuclear  devices  are 
known  to  have  been  used. 

The  numbers  in  the  columns  under  the  heading  'Processing 
Results'  refer  to  the  detection  status  of  the  events  as  recorded  at  each  station. 
A dash  indicates  that  the  event  was  not  processed  at  that  station--in  most 
cases  because  the  station  was  not  operational.  The  number  1 indicates  that 
the  event  was  detected,  2 indicates  that  it  was  not  detected,  3 indicates  that 
it  was  mixed  with  another  event,  4 indicates  no  data  were  recorded  for  that 
event,  and  5 indicates  that  a system  malfunction  obscured  the  data.  These 
codes  indicate  that  five  of  the  eighteen  events  listed  in  Table  VI- 1 cannot  be 
used  in  the  discrimination  study,  since  they  were  not  detected  at  any  of  the 
stations. 

Figures  VI-Z  to  VI-50  show  the  M^  - m^  plots  for  the  Seismic 
Research  Observatories  surface  wave  nnagnitudes  measured  at  periods  of 
15  seconds.  The  values  are  plotted  by  station  with  one  plot  encompassing  data 
from  all  regions  followed  by  individtial  plots  for  data  from  each  region.  Each 
regional  plot  also  shows  the  data  points  for  the  presumed  nuclear  explosions. 
The  symbols  used  in  these  plots  are: 

0 - earthquake  with  depth  less  than  60  km  or  depth  unknown 
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TABLE  VI- 1 

PRESUMED  NUCLEAR  EXPLOSION  EVENT  LIST 
WITH  DETECTION  STATUS 


I Detection  Statue* 

Event  Region  I ANMO 1 GUMO  I MAIO I NWAO  I TATO I SNZO 


* Meaning  of  Detection  Statue  Codee 

1 - event  detected  according  to  detection  critjeria  of 

Section  V 

2 - event  not  detected 

3 • mixed  event 

4 - eyatem  failure  reeulting  in  no  data  being  recorded 

5 - malfunction  (epike,  glitch  or  the  like)  in  aignal  gate 


ANMO  25.SECOND  M -nv  DATA  FOR  REGION  1 


ANMO  25-SECOND  M -m.  DATA  FOR  REGION 


GUMO  25-SECOND  M DATA  FOR  REGION  1 


GUMO  25-SECOND  M -nv  DATA  FOR  REGION  5 


GUMO  25-SECOND  . I -n^  DATA  FOR  REGION  7 


1W.I 


GUMO  25-SECOND  M -nv  DATA  FOR  REGION  12 


MAIO  2 5 -SECOND  M -m.  DATA  FOR  REGION  1 


^lO  25-SECOND  M -nrv  DATA  FOP  REGION  2 


MAIO  25-SECOND  M -m,  DATA  FOR  REGION  5 


MAIO  25-SECOND  M -m.  DATA  FOR  REGION  7 


MAIO  25-SECOND  M -m,  DATA  FOR  REGION  9 


SECOND  M -m,  DATA  FOR  REGION  1 1 


ici  ■■  ’iriii  J'l  ■ ■ ■ 


NWAO  25- SECOND  M -m^  DATA  FOR  ALL  RECttONS 
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NWAO  25-SECOND  M -m^  DATA  FOR  REGION  1 


NWAO  25 -SECOND  M -tn,  DATA  FOR  REGION  2 
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FIGURE  VI- 37 

TATO  25-SECOND  M -m  DATA  FOR  ALL  REGIONS 


TATO  25-SECOND  M -m,  DATA  FOR  REGION  1 


TATO  25-SECOND  M -m.  DATA  FOR  REGION  2 


FIGURE  VI-  41 

TATO  25-SECOND  M -nv  DATA  FOR  REGION  6 


TATO  25-SECOND  M -nx  DATA  FOR  REGION  7 


TATO  25-SECOND  M DATA  FOR  REGION  10 


TATO  25-SECOND  M -m.  DATA  FOR  REGION 


SNZO  25-SECOND  M -nx  DATA  FOR  ALL  REGIONS 


SNZO  25-SECOND  M -m.  DATA  FOR  REGION  5 


SNZO  25-SECOND  M -nv  DATA  FOR  REGION  12 


-f  - earthquake  with  depth  greater  than  60  km 

* - presumed  nuclear  explosion  from  Region  8 or  10 

- presumed  nuclear  explosion  from  Nevada  Test  Site. 

Such  depth  information  as  was  available  came  from  the  Nation- 
al Earthquake  Information  Service  bulletins.  The  straight  lines  in  each  plot 
represent  the  equation  of  the  M^-m^  relationship  for  that  data  set.  This  was 
computed  using  the  data  points  for  earthquakes  with  depth  less  than  60  km 
or  with  depth  unknown.  The  relationship  is  computed  with  a line-fitting  al- 
gorithm which  treats  neither  variable  as  dependent.  The  earthquakes  with 
depths  known  to  be  greater  than  60  km  were  excluded  from  the  fitting  pro- 
cedures to  avoid  biassing  the  fit.  (Deeper  events  tend  to  generate  relative- 
ly lower  values. ) The  lines  fitted  to  the  data  of  the  individual  regions 
are  based  for  the  most  part  on  rather  sparse  populations  and  are  intended 
to  be  used  only  as  a visual  aid  in  separating  earthquakes  from  presumed 
nuclear  explosions.  Little  emphasis  should  be  placed  on  their  particular 
slopes  and  intercepts. 

Referring  to  Figures  VI-2  to  VI-50  the  following  comments  can 
be  made  about  ^ose  events  of  Table  VI- 1 which  were  detected  at  one  or  more 
stations : 

a Event  13.  The  M values  for  this  event  as  recorded  at  ANMO 

s 

and  GUMO  show  excellent  separation  from  the  earthquake  pop- 
ulations, indicating  this  event  is  a nuclear  explosion. 

a Events  yl,  401.  410,  485,  500,  522,  523.  The  values 
for  these  Nevada  Test  Site  events,  while  in  general  on  the 
low  aide  of  ^a  aarthquake  population,  show  little  or  no 
separation  from  the  earthquake  population.  Therefore, 
if  the  only  means  of  differentiating  nuclear  explosions  from 


earthquakes  was  the  discriminant,  these  events  would 

be  classified  as  earthquakes. 


Event  537.  The  M values  for  this  event  as  recorded  at 
8 

ANMO,  GUMO,  MAIO,  and  NWAO  all  fall  in  the  earthquake 
population,  indicating  this  event  is  an  earthquake. 


Event  839.  The  M values  for  this  event  as  recorded  at 

8 

MAIO  fall  in  the  earthquake  population,  indicating  this  event 
is  an  earthquake. 


Event  1062.  The  M values  for  this  event  as  recorded  at 

s 

ANMO,  GUMO,  MAIO,  and  TATO  show  excellent  separation 
from  the  earthquake  population,  indicating  this  event  is  a 
nuclear  explosion. 


Event  1200.  The  M values  for  this  event  as  recorded  at 

8 

MAIO  and  TATO  show  excellent  separation  from  the  earth 
quake  population,  indicating  this  event  is  a nuclear  explo- 
sion. 


Event  1368.  The  transverse  component  M value  (the  only 

8 

component  on  which  it  was  detected)  for  this  event  as  recorded 
at  TATO  shows  excellent  separation  from  the  earthquake  pop- 
ulation, indicating  this  event  is  a nuclear  explosion. 


In  summary,  the  M^-m^  discriminant  appears  to  function 
quite  well  for  Eurasian  presumed  nuclear  explosions  but  not  for  Nevada  Test 
Site  presumed  nuclear  explosions.  This  discriminant  indicates  that  two  of 
the  detected  Eurasian  events  of  Table  VI- 1 are  earthquakes  and  the  remainder 
are  nuclear  explosions.  In  corroboration  of  this,  it  is  interesting  to  note 
that  almost  all  known  eastern  Kazakh  nuclear  explosions  have  m^  values 
greater  than  5. 0. 


D.  relationships 

Table  VI-2  presents  the  results  of  fitting  a straight  line 
to  the  M^-nri^  points.  The  M^-m^  relationships  are  presented  only  for  the 
case  of  all  data  points  measured  at  each  station.  The  M^-m^  relationships 
for  the  data  points  of  the  individual  regions  are  not  presented  since  in  most 
cases  the  number  of  points  is  too  low  to  yield  trustworthy  results. 

The  data  of  Table  VI-2  shows  that  only  at  NWAO  and  SNZO  is 
the  M^-m^  relationship  for  the  vertical  component  significantly  different 
from  the  M^-m^  relationship  for  the  transverse  component.  The  difference 
is  greater  at  SNZO.  where  fewer  M^-m^  points  were  available.  Therefore, 
the  M^-m^  relationships  at  these  two  stations  should  be  considered  only 
approximate  due  to  the  relatively  small  number  of  M^-m^  points  available. 


TABLE  VI-2 


SRO  STATIONS  M.-itib  RELATIONSHIPS 
FOR  25-SECOND  ENERGY 


n 


tation  Component 


AMMO 


i 


UMO 


MAIO 


NZO 


1.  65 

-3.83 

1.66 

-3.82 

1.41 

-2.88 

1.42 

-2.75 

1.33 

-2.21 

1.29 

-2.05 

-2.  93 
-3.36 


■2.  83 
•2.  91 


■3.78 

■5.08 


0.14 
0.  14 


0.  12 
0.12 


2 

4 


0.  09 

0.  11 


0.  13 
0. 18 


0.  06 
0.  07 


i I j| 


4.  92 
4.  90 


4.  87 
4.89 


4.53 

4.49 


5. 

5. 14 


4.55 
4.  54 


5.11 

5.39 


4.27  119 
4.3]  128 


156 

140 


3. 80  270 
3.77  268 


4. 66  60 

4.  64  59 


3.88  132 

4.15  141 


5.16  23 

4. 92  24 


Refer  to 
Figure  # 


VI- 2 


VI-11 


VI-19 


VI-29 


VI-37 


VI-46 


I where: 

M 

a 

= am^+b 

» 

i 

(t2 

= variance  normal  to  the  M^-m^  estimate 
and 

T 

1 

n 

= number  of  M points  used  in  compu^ 

tation. 

SECTION  vn 

SRO  NETWORK  DETECTION  CAPABILITY  ESTIMATES 


A.  DISCUSSION 


Up  to  this  point,  the  Seismic  Research  Observatory  stations 
have  been  evaluated  on  an  individual  station  basis.  This  section  investigates 
the  detection  capability  of  a network  composed  of  the  six  stations  under  eval- 
uation. Two  approaches  are  used  > the  first  utilizes  the  detection  statistics 
collected  for  the  individual  stations  and  the  second  utilizes  the  noise  level 
statistics  collected  for  the  individual  stations.  The  first  approach  is  used  to 
compute  detection  capability  estimates  for  the  short-period  network,  the 
long-period  network,  and  the  combined  short-period,  long-period  network. 
The  second  approach  is  used  to  compute  detection  capability  estimates  only 
for  the  long-period  network,  since  the  technique  used  does  not  have  the  cap- 
ability of  taking  into  account  the  effects  on  detection  capability  of  an  automatic 
detector. 


J 
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B.  NETWORK  DETECTION  CAPABIUTY  FROM  DETECTION  STATISTICS 
OF  INDIVIDUAL  STATIONS 


In  this  direct  estimate  of  network  detection  capability,  an 
event  was  declared  a detection  if  it  was  detected  at  one  or  more  stations  of 
the  network  using  the  detection  criteria  of  Section  V.  The  grouping  'ALL 
REGIONS  A'  and  'ALL  REGIONS  B ' have  the  same  meaning  as  in  Section  V . 
The  events  used  in  this  study  were  events  36  to  1384  of  the  data  base.  These 
events  were  selected  because  data  from  at  least  three  stations  had  been  pro- 
cessed for  each.  Just  as  in  Section  V,  the  50  percent  detection  threshold  is 
used  a-B  the  detection  capability  estimate. 
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The  detection  capability  estimates  for  three  networks  are 
presented  in  Figures  VII- 1 to  VII- 9.  First,  Figures  VII- 1 to  VII-3  present 
the  detection  capability  estimates  for  the  short-period  network.  Next, 
Figures  Vn-4  to  VII- 6 present  the  detection  capability  estimates  for  the 
long-period  network.  Finally,  Figures  VII-7  to  VII- 9 present  the  detection 
capability  estimates  for  the  combined  short-period,  long-period  network, 
where  an  event  is  considered  to  be  detected  if  it  is  detected  on  either  the 
short-period  or  long-period  network.  Each  of  the  'sub-figures'  in  Figures 
vn-1  to  Vn-9  represents  one  network  detection  capability  estimate.  The 
upper  portion  of  each  sub-figure  is  a histogram  giving  the  detection  statis- 
tics as  a function  of  bodywave  magnitude  for  that  particular  case.  The 
lower  portion  shows  the  percentage  of  events  detected  at  each  bodywave 
nnagnitude  (represented  by  asterisks),  the  fitted  maximum  likelihood  curve 
(represented  by  a solid  line),  and  the  90  percent  confidence  limits  for  this 
curve  (represented  by  dashed  lines).  The  values  for  'MB50'  and  'MB90' 
shown  on  the  figures  are  the  50  and  90  percent  detection  thresholds  as 
picked  from  the  maximum  likelihood  curve.  The  titles  for  each  sub-figure 
define  the  type  of  network  and  the  region  of  the  events  from  which  the 
detection  statistics  were  derived. 

Note  that  in  Figures  VII-l  to  Vn-9  the  terms  'ALL  REGIONS 
A'  and  'ALL  REGIONS  B'  differ  only  in  that  mixed  events  were  counted  as 
non- detections  for  the  detection  statistics  of  'ALL  REGIONS  A'  and  were 
rejected  from  the  detection  statistics  of  'ALL  REGIONS  B'.  (For  the 
events  used,  there  were  no  cases  in  vdiich  no  data  were  recorded  or  malfunc- 
tions obscured  the  data  at  all  operational  stations  simultaneously. ) Thus 
'ALL  REGIONS  A'  again  represents  the  absolute  detection  capability  estimate 
and  'ALL  REGIONS  B'  again  represents  the  conditional  detection  capability 
estimate  as  defined  in  Section  V.  For  short-period  data,  events  for  which 
the  automatic  detector  did  not  declare  a detection  are  counted  as  non-detec- 
tions. All  individual  region  detection  capability  estimates  are  conditional 
detection  capability  estimates. 
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FIGURE  Vn-3 

SHORT -PERIOD  NETWORK  DETECTION  STATISTICS 
REGIONS  10.  11,  AND  12 
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FIGURE  Vn-9 

COMBINED  NETWORK  DETECTION  STATISTICS 
REGIONS  10,  11,  AND  12 
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The  detection  capability  eatimates  of  the  three  networke  are 
summarized  in  Table  VII- 1.  The  following  points  should  be  noted  from  the 
data  in  this  table: 

e Short-period  detection  capability  is  poorer  than  long-period 

detection  capability  for  all  regions. 

e The  short- period  network  detection  capability  does  not 

suffer  from  mixed  events,  events  for  which  no  data  were 
recorded,  or  events  containing  malfunctions.  Thus,  the 
absolute  and  conditional  detection  capability  estimates  are 
the  same. 

# The  long-period  network  detection  capability  is  increased 

0.  2 m^  units  by  mixed  events.  Malfunctions  and  times  when 
no  data  were  recorded  at  individual  stations  had  no  effect  on 
the  long-period  network  detection  capability. 

e For  individual  regions,  the  short-period  and  long-period  net- 

works show  detection  capability  estimates  essentially  the 
same  as  the  lowest  estimate  for  that  region  from  the  individual 
station  estimates.  For  example,  the  long-period  network  esti- 
mate and  the  GUMO  station  estinaate  for  region  2 are  essentially 
the  same,  and  the  long- period  network  estimate  and  the  MAIO 
estlmata  for  region  11  are  essentially  the  same.  The  small 
differences  found  between  these  estimates  are  due  to  the  soirie- 
what  different  data  bases  used  to  generate  them. 

e The  contribution  of  the  short- period  network  to  the  combined 

network  is  to  replace  the  mixed  event  statistics  of  the  long- 
period  data  with  detected  and  non- detected  event  statistics 
of  the  short-period  data.  This  results  in  the  elimination  of 
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TABLE  VII- 1 

SRO  NETWORK  DETECTION  CAPABILITY  ESTIMATES 


1 50  Percent  Detection  Threshold 

Short-Period 

Net 

Long- Period 
Net 

Combined  Short- Period 
and  Long- Period  Net 

4.62 

4.  36 

4.  15 

4.  62 

4. 16 

4.  15 

4.93 

4.44 

4.51 

4.53 

4.  24 

4.20 

(4.2) 

4.  04 

3.  96 

4.31 

3.  98 

3.  96 

(4.4) 

3.  98 

3.89 

4.06 

3.86 

3.68 

(4.3) 

4.  10 

4.  00 

(4.7) 

4.  00 

4.05 

1 

(4.4) 

3.85 

3.97 

All  Regions  A 
All  Regions  B 
1 
2 

5 

6 
7 
9 

10 

11 

12 


All  Regions  A Mixed  events  were  counted  as  non-detections. 


All  Regions  B 
and  individual 
regions 

SP  Data 


Mixed  events  were  rejected  from  the  detection 
statistics. 


Events  for  which  the  automatic  detector  did  not 
declare  a detection  were  counted  as  non-detections 
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the  diiference  between  the  absolute  and  conditional  detection 
capability  estimates  observed  for  the  long-period  network. 

• For  individual  regions,  the  combined  network  shows  little 

difference  from  the  long- period  network.  The  differences 
observed  are  due  to  the  addition  to  the  detection  statistics 
of  events  which  were  mixed  on  the  long- period  data. 


C.  LONG- PERIOD  NETWORK  DETECTION  CAPABILITY  FROM  NOISE 
STATISTICS  OF  INDIVIDUAL  STATIONS 

This  approach  uses  a modified  version  of  the  network  capa- 
bility estimation  technique  reported  by  Snell  (Snell,  1976).  This  technique 
is  based  on  the  equations  presented  in  Table  VII- 2.  The  signal  amplitudes 
at  all  network  stations  are  calculated  for  an  assumed  event  at  epicenter  j 
with  magnitude  m^.  (equation  VII- 1).  Signal  and  noise  are  assizmed  to  be 
lognormally  distributed;  hence,  the  probability  that  the  signal-to-noise  ratio 
exceeds  the  station  detection  threshold  is  given  by  the  normal  cumulative 
probability  function  (equations  VII-2  and  VII- 3).  Signal  and  noise  are  meas- 
ured zero-to-peak  for  use  in  this  technique.  The  noise  mean  is  input  as  the 
geometric  mean  of  the  noise  amplitudes.  The  individual  station  detection 
probabilities  are  then  combined  into  the  network  probability  of  at  least 
station  detections,  0.(  > o ),  using  a recursive  relation  involving  the  individual 
station  detection  probabilities  (equation  VIl-4).  A detailed  description  of 
the  operation  of  this  technique  is  reported  by  Snell  (Snell,  1976). 

This  technique  was  modified  so  that  instead  of  computing  detec- 
tion threshold  magnitudes  (as  defined  by  Snell)  for  the  network,  a measure 
of  the  probability  of  detecting  an  event  of  given  magnitude  is  plotted.  This 
measure  of  the  probability  of  detecting  is  » where  P^  is 

the  probability  of  detecting  the  event.  This  value  was  chosen  as  the  measure 
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TABLE  Vn-2 

EVALUATION  OF  NETWORK  DETECTION  CAPABILITY 


1' 


M 

1 


I 

|] 

IJ 

II 


|] 

n 

El 


= m.  + t i y t E.. 

Pti 


* IO810SDT) 
2.2  1/2 


p.(>tt)  P.(k) 

J Icso  ^ 


(VU-1) 
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Symbols  above  are  defined  as  follows: 

• signal  amplitude  at  station  i for  event  j (0-p) 
- magnitude  of  event  j (m^) 


m. 


Mn 

2 

<^n 

4(x) 

N 


- standard  table  entries 

- station- epicenter  bias  corrections 

- mean  log^^  noise  amplitude  (0-p) 

- variance  of  log^^  noise 

- variance  of  ^oSjq  •ignal 

- normal  cumulative  probability  function 

- number  of  stations  in  the  network 


(k)  - probability  that  k stations  will  detect  event  j 

^^(>0)  - probability  that  0 or  more  stations  will  detect  event  j 

SDT  - station  detection  threshold,  i.  e. , signal-to-noise  ratio  re< 
quired  for  detection  at  station. 
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of  probability  of  detecting  rather  than  itself,  since  will  ro\ind  off  to 
one  for  high  probabilities  of  detection  (for  example  P^^  = 0.  999968  would 
be  rounded  off  to  1 > 0). 

The  mean  and  variance  of  the  ^<>82  q amplitudes  used 

to  produce  the  estimates  of  network  detection  capability  are  the  peak  25- sec- 
ond period  noise  amplitudes  described  in  Section  IV.  Station  reliability 
estimates  were  derived  from  the  station  reliability  estimates  and  probability 
of  mixed  events  estimates  of  Section  II.  The  signal-to- noise  ratio  required 
for  detection  at  a station  was  set  at  6 dB. 

Four  cases  of  long-period  network  detection  were  prepared 
for  this  report.  Figure  VII- 1 0 shows  the  probability  of  one  station  detect- 
ing an  5 event.  Figure  VII- 11  shows  the  probability  of  two  stations 

detecting  an  m^=4.  5 event.  Figure  VII- 12  shows  the  probability  of  one 
station  detecting  an  m^=5.  5 event.  Figure  VII- 13  shows  the  probability 
of  two  stations  detecting  an  mi^=5.  5 event.  Table  VII- 3 translates  the 
^10  ['  j^/l-Pp^  contours  into  probability  of  detection. 

It  is  difficult  to  make  comparisons  between  these  network 
detection  capability  estimates  and  the  detection  capability  estimates  pre- 
sented in  Table  VII- 1.  The  estimates  for  individual  regions  in  Table  VII- 1 
do  not  take  into  account  the  effects  of  mixed  events  while  the  estimates 
presented  in  Figures  VII- 10  to  VII- 13  do.  However,  some  qualitative  com- 
parisons can  be  made  by  asstiming  that  the  detection  thresholds  of  the  indi- 
vidual regions  will  b.e  raised  by  mixed  events  the  same  amount  as  were  the 
combined  regions.  This  would,  for  example,  make  the  long-period  net- 
work detection  capability  estimate  (50  percent  probability  of  detection)  for 
region  5 equal  to  4.  2 m^  from  Table  VU-l  and  the  90  percent  detection 
threshold  (from  Figure  VII-5)  equal  to  4.  5 m^.  Figure  VII- 10  shows  that 
the  probability  of  detecting  n Region  5 event  with  5 at  one  or  more  sta- 

tions is  approximately  0.  9. 
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m^  = 4.5 

2 Station  Detecting 


FIGURE  Vn-11 


PROBABILITY  OF  TWO  STATIONS  OF  THE  NETWORK 
DETECTING  AN  m.  =4.  5 EARTHQUAKE 
(LONg5>ERIOD) 
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FIGURE  Vn-12 

PROBABILITY  OF  ONE  STATION  OF  THE  NETWORK 
DETECTING  AN  m^«5.  5 EARTHQUAKE 

(LONG- PERIOD) 
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FIGURE  Vn-13 

PROBABILITY  OF  TWO  STATIONS  OF  THE  NETWORK 
DETECTING  AN  m=5.  5 EARTHQUAKE 
(LONG^ERIOD) 
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Thus,  there  appear*  to  be  fairly  good  agreement  between  the  two  methods 
of  estimating  network  detection  capability. 

Note  that  in  some  other  regions  the  indicated  standard  dev- 
iation of  the  event  population  is  very  large  in  contrast  to  the  assumption 
of  signal  standard  deviation  used  in  network  capability  evaluation  from 
noise.  These  may  be  due  to  large  variations  of  B-Factor,  source  radia- 
tion pattern  and  other  effects.  It  is  difficult  to  explain  large  variations  in 
observed  event  population  standard  deviations. 


SECTION  vin 
CONCLUSIONS 


This  •ection  summarizes  the  results  of  this  continued  eval- 
uation of  the  Seismic  Research  Observatory  stations.  The  major  results 
are: 

A.  DATA  QUALITY 

• In  general*  the  data  quality  was  good.  At  all  stations  except 
MAIO,  short-period  data  were  degraded  by  malfiuctions  one 
percent  of  the  time  or  less.  MAIO  short-period  data  was  de- 
graded by  a 1.7  Hz  waveform  of  unknown  origin  approximately 
seven  percent  of  the  time.  Long-period  data  were  degraded 
by  malfunctions  three  to  eight  percent  of  the  time. 

• Estimates  of  station  reliability  were  based  on  the  frequency 
of  occurrence  of  malfunctions  and  on  estimates  of  station 
down-time.  GUMO  had  the  lowest  reliability  (0.  7)  due  to 

a relatively  large  amount  of  station  down-time.  The  other 
stations  had  reliability  factors  of  0.  8 or  above. 

• The  probability  of  an  event  being  mixed  was  estimated  at 
approximately  0.  3,  with  little  variation  from  station  to 
station,  in  a range  of  m^  values  from  3. 1 to  5.  0. 

B.  SHORT-PERIOD  AUTOMATIC  DETECTOR 

• The  passband  filter  limits  placed  at  the  front-end  of  the 
short-period  automatic  detector  should  be  changed  from 
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0.  5-2.  0 Hz  to  0.  5-r<  5 Hz,  or  some  means  such  as  a notch 
filter  should  be  used  to  eliminate  the  1 . 7 Hz  waveform. 

• The  number  of  contiguous  points  for  which  the  short-term 
average  must  exceed  the  scaled  long-term  average  in  order 
for  a detection  to  be  declared  should  be  increased  from  20 
points  to  at  least  30  points  to  reduce  detection  declarations 
caused  by  glitches. 

• The  automatic  detector  appears  to  be  unstable  in  terms  of 
turning  off  after  a detected  event  dies  back  into  the  noise. 

C.  NOISE  ANALYSIS 

• In  order  of  increasing  magnitude,  the  mean  short-period 
RMS  noise  values  in  the  0.  5-4.  0 Hz  passband  (uncorrected 
for  instrument  response)  are:  0.  38mfx  at  ANMO,  0.  57  mji 
at  MAIO,  7.  69  m^  at  NWAO,  20.  61  m^  at  TATO,  28.  92  mfi 
at  SNZO,  and  40.  25  m^  at  GUMO.  These  values  show  a 
strong  correlation  with  station-nearest  coastline  separation 
indicating  that  a majority  of  short-period  noise  is  derived 
from  ocean  wave  energy  injected  at  the  coastline. 


• Based  on  mean  long-period  RMS  noise  levels  in  the  0.  023 

-0.  059  Hz  passband,  the  stations  can  be  divided  into  two 
groups--ANMO  and  MAIO  form  the  group  with  lower  RMS 
noise  levels  and  GUMO,  NWAO,  and  TATO  form  the  group 
with  higher  RMS  noise  levels.  Considering  the  locations  of 
these  stations,  this  suggests  that  there  is  at  least  some 
correlation  between  station- coastline  separation  and  long- 
period  noise  level. 
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The  long-period  RMS  noise  levels  at  SNZO  seem  too  high 
to  be  entirely  due  to  seismic  noise.  It  would  be  advisable 
to  check  the  instrximentation  at  this  station  to  determine 
whether  instrumentation  malfunctions  are  adding  to  the 
apparent  noise  levels. 

Oceanic  wave  energy  appears  to  play  a large  part  in  deter- 
mining the  noise  level,  spectral  content,  and  day-to-day 
variability  of  the  noise  field. 


DETECTION  CAPABILITY 

• No  short-period  events  of  the  data  base  were  detected  at 
SNZO,  while  only  two  were  detected  at  GUMO. 

• Mixed  events  had  very  little  effect  on  short-period  detection 
capability  estimates. 

• Station  down-time  caused  an  increase  of  approximately  0.  Z 
units  in  the  ANMO  short-period  detection  capability  estimates. 

• Malfunctions  caused  an  increase  of  approximately  0. 1 

units  in  the  MAIO  short-period  detection  capability  estimates. 

• Mixed  events  caused  an  increase  of  approximately  0.  3 m^ 
units  in  the  NWAO  long-period  detection  capability  estimates. 

• Mixed  events,  system  failures,  and  malfunctions  equally  con- 
tributed to  an  increase  of  approximately  0.  3 m^  units  in  the 
SNZO  long-period  detection  capability  estimates. 

• Mixed  events,  and  to  a lesser  degree  system  failures  and 
malfunctions,  caused  an  increase  of  approximately  0. 4 
units  in  the  long-period  detection  capability  estimates  of 
ANMO,  MAIO,  and  TATO. 
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• System  failures,  and  to  a lesser  degree  mixed  events, 

caused  an  increase  of  approximately  0.  5 m units  in  the 

b 

GUMO  long-period  detection  capability  estimates. 

E.  DISCRIMINATION 

e Surface  wave  magnitudes  of  presumed  nuclear  explosions 

from  the  Nevada  Test  Site  fall  in  the  earthquake  population 
of  the  M^-m^  discriminant. 

• Events  537  and  839  show  surface  wave  magnitudes  indicating 
that  they  are  earthquakes. 

• Events  13,  1062,  1200,  and  1368  show  surface  wave  magnitudes 
indicating  that  they  are  nuclear  explosions. 

F.  NETWORK  DETECTION  CAPABILITY 

• Using  direct  estimates  of  detection  capability  derived  from 

the  events  of  the  data  base,  it  was  found  that  the  combined  short- 
period,  long-period  network  shows  little  difference  from  the 
long-period  network  detection  capability  estimates.  The 
primary  contribution  of  the  short-period  network  is  to  replace 
the  mixed- event  statistics  of  the  long-period  data  with  detect- 
ed and  non-detected  event  statistics  of  the  short-period  data. 

• Short-period  detection  capability  is  poorer  than  long-period 
detection  capability  for  all  regions.  This  appears  to  be  due 
to  the  effect  of  automatic  detector  on  the  short-period  data. 

• The  short-period  network  detection  capability  estimates  are 
not  affected  by  mixed  events,  malfunctions,  or  system  failures 
at  individual  stations. 
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• The  long'period  network  detection  capability  estimates 
are  increased  0.  2 units  by  mixed  events.  Malfunctions 
and  system  failures  at  individual  stations  did  not  affect  these 
estimates. 

• For  individual  regions,  the  short-period  and  long-period 
networks  show  detection  capability  estimates  essentially  the 
same  as  the  lowest  estimate  for  a given  region  from  the  in- 
dividual station  estimates. 

• Based  on  such  comp>arisons  as  are  possible,  the  above- 
described  long-period  direct  network  detection  capability 
estimates  appear  to  agree  quite  well  with  the  network  detec- 
tion capability  estimates  derived  from  noise  statistics. 

Future  work  on  the  evaluation  of  the  Seismic  Research  Ob- 
servatory stations  should  be  directed  toward  the  following  points: 

• Evaluate  all  stations  as  they  become  operational. 

• Expand  the  short-period  and  long-period  noise  data  bases  so 
that  at  least  one  year  of  data  is  available  at  each  station  for 
the  study  of  noise  level  trends. 

• Obtain  continuous  short-period  data  and  the  automatic  de- 
tector algorithm  so  that  the  effects  of  changes  in  the  detec- 
tor input  parameters  can  be  studied. 

9 Expand  the  signal  data  base  so  that  detection  and  discrimin- 

ation results  can  be  determined  for  new  stations  and  better 
defined  for  particular  regions  for  stations  already  under 
evaluation.  This  would  also  be  used  to  better  define  station 
reliability  and  mixed  event  probability  estimates. 

• Investigate  the  discrimination  of  earthquakes  and  presumed 
nuclear  explosions  using  short-period  data. 
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APPENDIX  A 
THE  DATA  BASE 


The  data  bate  ueed  for  this  evalxiation  of  the  Seismic  Research 
Observatory  stations  is  presented  on  the  following  pages.  Table  A-1  lists 
the  events  from  the  data  base  used  in  the  evaluation  of  short-period  and  long- 
period  data  recorded  at  each  station. 

In  the  data  base  list,  the  column  headed  'EVNO'  gives  the 
identifying  number  assigned  to  each  event.  The  columns  headed  'DATE'  and 
'TIME'  give  the  date  and  origin  time  of  each  event.  The  epicentral  coordin- 
ates are  listed  under  'LAT.  ' and  'LONG'  for  latitude  in  degrees  north  and 
longitude  in  degrees  east.  The  bodywave  magnitude  is  listed  under  the  head- 
ing 'MB'.  The  column  headed  'Q'  gives  the  NORSAR  quality  rating  for  the 
parameterit!  of  each  event  where 

1 3 good  to  excellent  quality 

2 3 fair  to  good  quality 

3 3 poor  to  fair  quality. 

The  seismic  source  area  is  listed  in  the  column  headed 
'LOCATION'.  Finally,  the  seismic  sub- region  as  defined  by  Flinn  and 
Engdahl  (Flinn  and  Engdahl,  1965)  for  each  event  is  listed  under  the  heading 
'SUBREG'. 
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Station 

1 Event  Nos.  of  Processed  Events) 

SP 
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1-1402 
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APPENDIX  B 

INSTRUMENT- RESPONSE  CORRECTED  NOISE  CHARACTERISTICS 


t,-. 


1^ 


In  the  course  of  the  evaluation  of  the  Seismic  Research  Obser- 
vatory stations,  instrument- response  corrected  noise  field  characteristics 
were  studied  at  each  station  under  evaluation.  In  this  appendix  to  the  con- 
tinuation of  the  Seismic  Research  Observatories  evaluation,  the  extension 
to  this  study  is  presented. 

To  correct  the  data  for  the  effects  of  the  instrument  response, 
the  data  were  first  Fourier  transformed  to  the  frequency  domain.  Noise 
amplitude  spectra  were  then  computed  and  smoothed  to  64  frequencies.  The 
smoothing  was  performed  to  simplify  the  task  of  making  the  instrument  re- 
sponse corrections.  In  the  course  of  this  smoothing,  the  spectral  values 
were  scaled  so  that  Parseval's  formula  would  hold.  This  formula,  which 
permits  the  computation  of  RMS  noise  from  frequency  domain  data,  is  ex- 
pressed as: 


RMS'*  = 

a 


|A(f.)|^*C(f.)‘ 
i=a  ^ 


1/2 


(frequency-domain 

estimate) 


where 


|A(f.)l 


C(f^) 


a 

b 


- the  elemental  frequency  interval, 

= the  discrete  Fourier  transform  spectral  density  es- 
timate at  frequency  f^, 

= the  instrument  response  correction  at  frequency  f^, 

- the  initial  frequency  index 
s the  final  frequency  index. 
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I ‘ei 


During  the  course  of  studying  the  noise  field  this  year,  the 
RMS  noise  values  produced  by  this  method  were  checked  against  those  pro 
duced  by  bandpass  filtering  the  data  and  computing  RMS  noise  from  the 
equation  i — — • 


RMS 


(time -domain  estimate) 


number  of  data  points,  and 
.th  . 

the  1 data  point. 


Since  this  equation  does  not  allow  for  making  instrument  re- 
sponse corrections,  C(f  J in  Parseval's  formula  above  was  set  equal  to  one. 
The  RMS  values  so  computed  with  Parseval's  formula  were  found  to  be 
approximately  10  percent  higher  than  those  computed  with  the  second  equa- 
tion. Therefore,  before  proceeding  further,  it  was  necessary  to  investigate 
the  causes  of  this  difference  between  the  frequency-domain  and  time-domain 
estimates  of  RMS  noise. 


The  algorithm  used  to  filter  the  data  was  the  first  point  of  in- 
quiry. The  design  of  this  filter  is  such  that  if  one  specifies  a passband  of 
0.  023-0.  059  Hz  with  0.  002  taper,  the  corners  are  not  at  0.  023  Hz  and  0.  059  Hz 
as  might  be  expected  but  are  0.  02466  Hz  and  0.  05737  Hz.  (This  is  illustrated 
by  Fibure  B-1.  ) Thus,  instead  of  passing  all  energy  with  periods  between  16.  9 
and  43.  5 seconds  as  intended,  the  filter  passes  all  energy  with  periods  between 
17.4  and  40.  6 seconds.  The  difference  at  the  shorter-period  end  is  insignifi- 
cant. At  the  high  end,  however,  this  difference  results  in  more  energy  being 
added  to  the  RMS  estimate  than  desired. 


The  next  stage  in  investigating  the  difference  between  the  fre- 
quency-domain and  time-domain  RMS  estimates  was  to  test  the  two  methods 


on  synthetic  data  consisting  of  sine  waves  of  varying  periods.  This  data  was 
processed  using  periods  between  10  and  50  seconds,  with  1 second  increments. 
The  results  shown  in  Table  B-i  in  the  column  headed  'T.  D.  ESTIMATES'  are 
exactly  what  would  be  expected--  outside  the  passband  limits,  the  small  non- 
zero values  are  due  to  aide- lobes  produced  by  the  Fourier  transform.  The 
values  inside  the  passband  limits  are  simply  the  RMS  value  of  a sine  wave. 

The  increasing  values  between  the  periods  of  16  and  18  seconds  and  the  de- 
creasing values  between  the  periods  of  40  and  44  seconds  are  due  to  the 
cosine- squared  taper  of  the  filter.  ■ 

The  results  of  applying  Parseval's  formula  without  instrument 
response  corrections  to  the  synthetic  data  after  Fourier  transforming  and 
smoothing  to  64  frequencies  is  shown  in  Table  B-1  under  the  heading  '64 
FREQ  SMOOTHING- -F.  D.  ESTIMATE'.  The  results  here  are  the  same 
as  the  time-domain  estimates  at  the  shorter  periods.  However,  at  longer 
periods,  the  frequency-domain  estimates  do  not  drop  off  until  47  seconds 
period,  indicating  that  too  much  energy  is  being  added  into  the  RMS  noise 
estimate. 

This  difference  at  longer  periods  is  due  to  the  smoothing  pro- 
cedure. By  smoothing  the  spectra  to  64  frequencies,  the  energy  is  divided 
into  bins  of  width  0.  0039062  Hz.  These  bins  are  too  wide  to  permit  the  fre- 
quency domain  RMS  values  being  the  same  as  the  time  domain  values.  After 
some  trial-and-error  tests,  it  was  found  that  if  the  spectra  were  smoothed 
to  128  frequencies  instead  of  64  frequencies,  the  frequency-domain  RMS 
values  would  differ  from  the  equivalent  time-domain  RMS'Values  by  less 
than  2 percent.  This  difference  was  considered  to  be  acceptable,  so  instru- 
ment-response corrected  RMS  noise  values  were  computed  using  Parseval's 
formula  and  noise  spectra  smoothed  to  128  frequencies. 

Figures  B-2  through  B-7  show  the  instrument- response  cor- 
rected RMS  noise  values  plotted  against  Jxilian  day,  udiile  Figures  B-8  through 
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TABLE  B-1 
RMS  NOISE  TEST 

BP=0.  023-0.  059  WITH  0.  002  TAPER 
I/P=SINE  WAVE  WITH  MAX  AMP=1.0 
NO  INSTRUMENT  RESPONSE 


64  Frequency 
Smoothing 


r requency 
Smoothing 


riod 
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INSTRUMENT -RESPONSE  CORRECTED  RMS  NOISE  VALUES 
PLOTTED  AGAINST  JULIAN  DAY  FOR  MAIO 
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INSTRUMENT -RESPONSE  CORRECTED  RMS  NOISE  VALUES 
PLOTTED  AGAINST  JULIAN  DAY  FOR  SNZO 


B-13  show  the  RMS  noise  level  trends  computed  from  these  values.  Since 
these  figures  differ  from  those  of  Section  IV,  where  the  instrument  response 
was  not  corrected  for,  only  in  level,  no  further  discussion  will  be  presented. 

The  mean  values  and  associated  standard  deviations  of  the  RMS  noise  are 
presented  in  Table  B-2. 

Figures  B-14  through  B-19  present  the  mean  instrument- re- 
sponse corrected  RMS  amplitude  spectra  for  each  of  the  six  stations  under 
evaluation.  The  left  half  of  each  figure  shows  the  log  of  the  RMS  amplitude 
spectra  with  vertical  bars  showing  plus-or-minus  one  standard  deviation. 

The  mean  noise  spectra  presented  for  ANMO  in  the  preliminary 
evaluation  showed  much  higher  noise  levels  on  the  north  component  than  on 
the  other  two  components  for  periods  above  30  seconds.  This  is  not  the  case 
in  Figure  B-14.  The  difference  is  due  to  the  rejection  from  the  noise  data 
base  this  year  of  noise  samples  containing  glitches  on  the  north  component,  some 
of  which  were  included  in  the  noise  data  base  last  year. 

The  mean  noise  spectra  presented  last  year  for  noise  recorded 
at  GUMO  displayed  a local  maximum  at  26  seconds  period,  a feature  which 
does  not  appear  in  Figure  B-15.  The  disappearance  of  this  feature  is  probably 
due  to  the  averaging  together  of  the  spectra  of  many  more  noise  samples  this 
year.  Whatever  the  cause  of  this  feature,  it  apparently  died  out  in  the  months 
following  April,  1976.  The  remaining  mean  noise  spectra  shown  in  Figures 
B-16  through  B-19  are  not  significantly  different  from  those  presented  in  the 
preliminary  evaluation. 
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FIGURE  B-8 

RMS  NOISE  LEVEL  TRENDS 
PLOTTED  AGAINST  JULIAN  DAY  FOR  ANMO 
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FIGURE  B>12 

RMS  NOISE  LEVEL  TRENDS 
PLOTTED  AGAINST  JULIAN  DAY  FOR  TATO 
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FIGURE  B-13 

RMS  NOISE  LEVEL  TRENDS 
PLOTTED  AGAINST  JULIAN  DAY  FOR  SNZO 


TABLE  B-2 

MEAN  RMS  INSTRUMENT-RESPONSE  CORRECTED  NOISE 


Station 


ANMO 

GUMO 

MAIO 

NWAO 

TATO 

SNZO 


11.71  4.72 

13.27  4.  13 

9.  72  5.  35 

16.63  6.00 

16.93  7.89 

49.33  16.56 


Mean  S. 


10.64  6.03 

20. 78  6.  84 

9.19  4.19 

21.31  7.28 

18.44  9.13 

30.23  12.66 


Mean 

S.D. 

No.  Samples 

11.87 

4.  97 

221 

21.44 

6.  32 

151 

9.  58 

5.01 

207 

14.56 

6.25 

167 

21.32 

11.84 

95 

34.56 

14.  52 

141 

SITE  36 


